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Abstract 



This paper is a continuation of [Brlvr] and [Ivrl, Ivr5, Ivr6]. I derive sharp spectral 
asymptotics (with the remainder estimate 0([i~ 1 h 1 ~ d + /U 2 _1 /i 1_ 2 ) for <i-dimensional 
Schrodinger operator with a strong magnetic field; here h and \i are Plank and binding 
constants respectively and magnetic intensity matrix has full rank at each point. 

Introduction 

0.0 Preface 

In this paper I consider multidimensional Schrodinger operator 
(0.1) A = A + V(x), A = V P 3 g' k (x)P k , 



j,k<d 



p. = hDj - nVj{x), h e (0, 1], n > 1. 
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It is characterized by magnetic field intensity matrices (Fjk) with Fjk = (d x ,Vk — d Xk Vj), 
which is skew-symmetric d x ci-matrix, and (F£) = (g^ k )(Fk p ) which is equivalent to skew- 
symmetric matrix (g^ k )^ (Fjk)(g^ k )^ . Then all eigenvalues of (F£) are idf p (fi > 0, p = 
1, . . . , r) and of multiplicity d — 2r where 2r = rank(i^). 

In this paper I assume that magnetic field intensity matrix has a full rank: 

(0.2) rank(i^) = 2r = d, < c o 

and (under certain conditions) I derive sharp spectral asymptotics (with the remainder 
estimate 0(n~~ 1 h l ~ d + /i2 -1 /i 1- 2 J). The typical example (already studied in [Ivrl, Ivr6]) is 
2-dimensional magnetic Schrodinger operator. 

I consider operator in some domain or on some manifold X with some boundary condi- 
tions, assuming that it is self-adjoint in L 2 (X) and denote by e(x,y,r) Schwartz' kernel of 
it spectral projector. 

0.1 Canonical form 

It is known (see e.g. [Ivrl], Chapter 6) that if X = M. d and g^ k ,Fjk,V are constant then 
operator is unitary equivalent to 

(0.3) E f j (hW j + Sx}) + E h2D ' + v 

l<j<r r+l<j<d-2r 

and e(x,x,r) = ^ w (r) where under condition (0.2) ^ 

(0.4) £^(r) ^ (2vr)-V^- r E °{ T ~ E^ + Wi** ~ V )h ■ • ■ fry/g 

with g = det(gjk), (gjk) = (fi^) -1 and Heaviside function 9. 

In particular, under condition (0.2) (and only in this case) spectrum is pure-point (of 

infinite multiplicity) consisting of Landau levels E a = ^2j(2ocj + l)fj, a E Z +r . 

Now in the general case (for variable V and, may be, g^ k , Fjk) I am interested in asymp- 
totics of f if)(x)e(x, x, 0) dx as h — > +0, ji — > +oo where tf> is a fixed function, smooth and 

^For 2r < d 

(0.4)* £ d MW (r) d ^c d _ 2l .(2^)" d +>^- rf +'' E ( r - E( 2a J + X )f^h - V^'h ...f ry /g 

where u>k is a volume of unit ball in R fc and 2r = rank(i 7 jj ; ). 
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compactly supported in X. The principal part of such asymptotics under some reasonable 
and obvious conditions is of magnitude h~ 2r + fi r h~ r . On the other hand, the remainder es- 
timate is 0(fih l ~ 2r + // h~ r ) unless one imposes some non- degeneracy condition while under 
the strongest possible condition one can expect remainder estimate 0(/i _1 /i 1_2r +/i r ~ 1 /i 1 ~' r ). 

It is known [Ivrl] that in the smooth case for d = 2 one can reduce operator in question 
to canonical form 

(0.5) Yl ^~ 2k ~ 2m ~^rn,kA^^- 1 hD 2 )(h 2 Dl + fx 2 xl) m h^ 

m+k>l 

with dip = F o v|/ , q 01 = y o \J/ and certain diffeomorphism \l>o : TM} — > IR 2 . Then one 
can replace harmonic oscillator h 2 D\ + [i 2 x\ by one of its eigenvalues (2a + l)/ih; we ignore 
terms with j > 1. 

So, for \ih < e our operator looks like a family (with Co(/i/i) _1 elements) of 1-dimensional 
/i~ 1 /i-pdos with principal symbols ((2a + l)jihF + V) o \[> . Spectral asymptotics for 
them is very sensitive to a degeneracy of the symbol but under non-degeneracy condi- 
tion |6| + |V6| > ei, b = ao,i + ai,o it has a principal part of magnitude nh~ l and a 
remainder 0(1). This non- degeneracy condition is equivalent to | VV/F\ > e± and the final 
asymptotics has its principal part of magnitude h~ 2 and a remainder 0()Lt _1 /i _1 ). 

As fih > €q we have the family of no more than Cq of /i~ 1 /i-pdos, non- degeneracy 
condition transforms into \(2at + l)fih{V/ F) \ + |V(V/i^)| > e\ and the final asymptotics 
has its principal part of magnitude /i/i -1 and a remainder 0(1). 

Non-smoothness prevents us from complete canonical form but we have "poor man" 
precanonical form which is sufficient [Ivr6]. 

Multidimensional case is much more tricky. Under condition (0.2) one could expect 
canonical form 

(0.6) ^ k - 2 ^a m>k {x'\ ^hD")H m 

m£l+ r ,k:\m\+k>l 

with H = (Hi, . . . ,H r ), harmonic oscillators Hj = h 2 D 2 + /i 2 x| and a( mj o) = fj ° as 
m k = S jk , a^) = V o^ , diffeomorphism # : TW = R 2r -> R 2r , x" = {x r+u x 2r ). 

In this case for fih < e one would get a family (with (fih)~ r elements) of r-dimensional 
/i^/i-pdos and under proper non-degeneracy condition |V(/j/V)| > e-y the principal part 
of asymptotics for them would be of magnitude [i r h~ r and a remainder 0(fi r ~ 1 h 1 ~ r ) and 
the final asymptotics has its principal part of magnitude h~ 2r and a remainder 0(/i _1 /i 1_r ). 

As \xh > €q one would have the family of no more than Co of /^ _1 /i-pdos, under proper 
non- degeneracy condition \b\ + V6| > e\, b = + ty^hfj + V, the final asymptotics 

has its principal part of magnitude fi r h~ r and a remainder 0(n~ 1 h 1 ~ r ). 
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0.2 Resonances 

However, the resonances prevent us from reducing our operator to the desired canonical 
form. In fact, m-th order resonances 

(0.7) hli = as 7 6 Z r and M d = ^ W = m > 2 

3 3 

prevent us from reducing properly terms 

(0.8) Yl ^- k - H - m a a , k {x\^- l hD l ')(hD') a {^'f. 

a,P m .\a\ + \p\=m; k 

In particular, due to the second order resonances we can reduce the main part of operator 
only to 

(0.9) ^ ajk(x", n~ l hD")Z*Z k + a (x", uT l hD") 

mean j,fcem 

with Zk = hDk + ifiXk where m G 971 are disjoint subsets of {1, . . . , r} and eigenvalues 
of each of matrices (ajk)j,kem are close to one another. This leads to the necessity of the 
matrix rather than the scalar non- degeneracy (microhyperbolicity) condition. 

Further, the third order resonances fi = fj+fk (with possible j = k) prevent us from get- 
ting rid off the terms (f~ 1 a,jf lt t(x", \i~ x hD")Z*Z k Zi and their adjoints ar~ x a*- kl {x" , n~~ x hD")Z 'jZ£Z* 
and these terms appear as perturbations of the main part of operator unless Fj k and g^ k 
are constant. 

The fourth order resonances fi = fj + fk + fi and + fj = fk + fi leave us with the 
whole bunch of terms instead of /i _2 if| but these terms are smaller than those produced 
by the third order resonances and so on. 

More precisely, the symbol of the reduced operator is 

(0. 10) n 2 Y Yl a JkCjCk + a + [i 2 Y Yl Re a M(j(k(i + ■■■ 

meOTj,fcem ne^tj,k,len 

where Q = £j + iXj, ' means complex conjugation for scalars and Hermitian conjugation 
for matrices; precise definition of 9t will appear later. 

In order to recover the desired remainder estimate we need to consider propagator trace 

Tr(e _ih 1<j4 ) for \t\ < T\ = f e/i with small constant e. 

Note that the Poisson brackets of this symbol with a m = f ^ . €m a^fcCjCfc are of magnitude 
1 and therefore variation of a m at time interval [— T 1; 7\] is not necessarily small. 
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On the other hand, higher order terms have 0(fi 1 ) gradients and therefore the Poisson 
brackets of (0.10) with 

(0.11) a° n = 5> ifc CjC* + ^ _1 E Re %«C]CfcC/ 

are of magnitude and therefore variation of a n at time interval [— T\,Ti] is small. 
Therefore, in the microhyperbolicity condition we need to consider a n rather than a m ; 
however, the third-order terms actually do not appear in this condition. 

0.3 Microhyperbolicity 

Thus we arrive to the final microhyperbolicity condition at point z = (a/,f ') and thus at 
point x = ( z )- 

Condition 0.1. Assume that 

(0.12) First of all we assume that set {l,...,r} is split into disjoint groups m G Wl(z) 
and matrix (a^) is block-diagonal and in a small vicinity of z its blocks are close to scalar 
matrices / m i# m where / m are disjoint. 

(0.13) Further, we assume that set {1, . . . ,r} is split into larger disjoint groups n G %l(z) 
such that no intergroup third order resonance is possible: 
\fi ~ fj ~ fk\ < eo i, j, k G n for some n. 

(0.14) Finally, we assume that for each r = (r n ) ne «n such that | Xlne^n 7 "" + ^1 — e there 
exists vector I = £(z, r) G M. 2r such that 

(0.15) /i 2 ^(a;/iV))CjG>ei 
as long as 

(0.16) a i*CjCfc = An VC G C r Vn G 9t 

To impose this condition we need to assume first that a does not vanish: 
(0.17) |V| > e . 

Condition 0.2. This microhyperbolicity condition should be modified as ah > e^: in 
(0.15) — (0.16) one needs to consider only r n = ^ ien (2ai + l)fihfi with ccj G Z + ; now r 
runs a discrete set. 
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On the other hand, for constant condition (0.15)-(0.16) is equivalent to 
(0.18) |W|>e x ; 
and for fih > one should modify it to 

(0.19) \'^2(2a i + l)fih + V\ + \VV\>e 1 Va G Z +r . 

i 

Then we basically do not need condition (0.17). Also in this case £ does not depend on r. 

Situation becomes even more tricky due to non-smoothness of coefficients. Let us list 
our assumptions first. 

0.4 Assumptions and results 
0.4.1 Assumptions 

We assume that in -8(0, 1) the following conditions together with (0.2), (0.17) hold: 
(0.20) c-'Kj^Q^k^c V£:|£| = l, 

jk 

(o.2i) /eC^(s(o,i)), 

(0.22) Vj = drfj, (f>j G C l+2 ~ a 

(0.23) V G C l ' a 

where C l,a for is a class of functions such that its Z -derivatives are continuous with conti- 
nuity modulus t l ~ l °\ logt| _<J where Iq = [l\ unless / G Z + , o < in which case — 1 — 

I suspect that condition (0.22) is rather technical and one can replace it by a more 
natural assumption 

(0.22)' F jk G G l >° 

Further, in the non-smooth case our theorems deal not with spectral kernel e(x, y, r) of 
the original operator A but with spectral kernel e(x, y, r) of it framing approximations A 

2 )One can take two copies of approximation such that our estimates hold for both and A~ < A < A + 
in the operator sense. We call them framing approximations. 
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0.4.2 Results: weak magnetic field 

The first statement deals with a (respectively) weak magnetic field: 

Theorem 0.3. Let conditions (0. 2), (0. 11), (0.20) -{0.23) with (I, a) h (I, a) h (1,2) and 
condition 0. 1 be fulfilled and 

(0.24) v>eCp(£(o,~)), o<v<i. 

Let /i < h 5 ^ 1 with arbitrarily small exponent 5 > 0. 
(%) Assume that one of the following two conditions is fulfilled: 

(a) (/,*)>: (2,0) 3 > 

(b) #W= 1 and (I, a) >z (1,2). 
Then the following estimate holds: 

(0.25) U=\ J (e(x,x,0) -S MW (x,0)^(x)dx\ < 

Cii~ x ti- d + C(ph) l \ logh\ l - a h- d + Cfi{fih) l+1 \ \ogh\ l+1 -~ a h- d ; 

(ii) In the general case for (/, a) -< (2, 0) the following estimate holds: 

(0.26) TZ < Cii l - [ h l ~ d \ log/i|- CT + C(ph) l \ logh\ l - a h- d + Cfi{fih) 1+1 \ log h\ l+1 -~ a h~ d . 

Corollary 0.4. (i) For smooth operators the best remainder estimate 

(0.27) TZ < c^h 1 -* 

holds as ji < h 6 ' 1 with arbitrarily small exponent 5 > 0. 

(ii) If (I, a) >z (3,1) then asymptotics (0.27) holds for ji < Ch~2 | log h\~2 . 

(Hi) For fi < Ch~z \ \ogh\~z one can skip last term in estimates (0.25) — (0.26). 

3 \l, a) y (V, a') if either I > I' or I = I', a > a' etc. 
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0.4.3 Results: intermediate magnetic field 

Our second statement deals with a stronger magnetic field: 

Theorem 0.5. Let conditions {0.2), (0.17), {0.20) -{0.23) with (J, a) y {I, a) h (1,2), (l,a) h 
(2,1) and condition 0.1 and (0.24) be fulfilled. 
Then for 

(0.28) nl = Chr*\]ogh\-* < n < fi* 2 = e'/r 1 ! log/i]" 1 

with sufficiently small constant e' > the following estimate holds: 

(0.29) Tlx = | J (e(x,x,0) -S MW (x,0) - E^ r (x,0)^(x)dx\ < 

C^h 1 " 1 + Cfi~ l \ log/i|~ CT /r d 

where formula for non-smoothness correction Sf^ r is rather complicated and given in sub- 
section 5.2.4 but one can take it unless I > 2, (I, a) -< (3,0) when it does not exceed 
Cfi~ 2 (fih) l ~ 2 h~ d \ log(/i/i)|- CT (for l>2). 

Remark 0.6. (i) Condition (0.28) means exactly that estimate (0.29) is stronger than 
(0.25), (0.26). 

(ii) If (I, a) y (3, 1) then the best remainder estimate (0.27) holds as fi < \i* 2 . 
0.4.4 Results: strong magnetic field 

Theorem 0.7. Let conditions (0.2), (0.17), (0.20) -(0.23) with (l,a) y (I, a), (I, a) y (2,1) 
and Condition 0.1 and (0.24) be fulfilled. 
Let 

(0.30) V* 2 <Cn< fi* 3 = e'h" 1 

with sufficiently small constant e' > 0. Then 

(i) //#01 = 1 then the following estimate holds: 

(0.31) Tlx < Cfi^h 1 ^ + C(fi~ l h)^\\ogh\^ a h- d ; 

in particular, for (I, a) y (2,1) remainder estimate 0(ix~ 1 h l ~ d ) holds. The description of 
^lcorr i s same as i n 0.5. 

(ii) If I > 2 then estimate (0.27) holds. 

(Hi) If > 1, I < 2, then the following estimate holds: 

(0.32) Tl x < Cfi^h 1 ^ + Cfi~ l \ logh\- a h- d e Collhlloehl . 
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Remark 0.8. (i) For #91 = 1 estimate (0.31) remains valid under relaxed restriction to /x: 

(0.30)' /i* < C/i< /4* = C7T 1 . 

For yU > \x\ one can replace Condition 0.1 by Condition 0.2. 

(ii) If F 3 k = const then #91 = 1, microhyperbolicity condition is equivalent to (0.18) and 
one can skip condition (0.17). 

0.4.5 Results: superstrong magnetic field 

Now we need to consider case 

(0.33) [i > /4* = e'h" 1 

in which spectral gaps are possible. I will discuss them in the end. 

In this case we need to modify our conditions. First of all one must assume that #91 = 1 
which will be covered by other conditions I impose. 

Second, we need to replace (0.23) by 

(0.34) 3a = aeZ +r : V + ^(2a 3 - + l)fihfj G C l '° '. 

3 

Further we need to assume that one of the following conditions is fulfilled: 
(0.35) g jk = const, F jk = const and (0.19) holds 
(0.36) a in (0.34) is unique (at each point): 

(0.37) Va G Z +r , a ^ a | ^(^ - a^/il > e o > 0; 

3 

in the former case we impose microhyperbolicity condition (0.19); in the latter case we 
impose the following microhyperbolicity condition: 

(0.38) |v(v(^(2a, + l)/ J )" 1 )|>6 . 

3 

Theorem 0.9. Let conditions {0.2), (0.20), (0.21) -(0.22), {0.S '4) with (I, a) y {l,<r), h 
(2,1) and (0.24) be fulfilled. 
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Finally, let either conditions (0.36), (0.19) or (0.37), (0.38) be fulfilled. Then 

(0.39) K 2 =| J(e(x,x,0)-S MW (x,0)-S™Z(x^)^(x)dx\ < 

Cii r - l h}~ r + C/j J r h~ r (ii~ 1 h)i\\ogfi\i-' J + C// r+1 /i 1 - r (/i~ 1 /i)5| l g/i|5-^ 

with 

(0.40) ^(x,r) d = f (27r)->V-'' £ ^-^(2^ + 1)/^-^' 

wii/i H 7 ' = 0(1) (the exact formula is rather complicated). 

Finally, consider spectral gaps: 

Theorem 0.10. Let conditions (0.2), (0.20), (0.21) - (0.22), (0.34) with (I, a) h (I, a), 
(I, a) y (2,1) and (0.24) be fulfilled. 

Further, let either condition (0.36) be fulfilled or \ih < C\. Finally, let 

(0.41) | J2( 2a i + l )» h + ^1 > e i Va e Z+r - 

i 

Then 

(0.42) K < Cfi~ s 

with arbitrarily large s. 

Surely, microhyperbolic case of theorem 0.9 is separated from elliptic case of theorem 
0.10 by degenerate one but we do not consider it (it is possible assuming that the degeneracy 
is not (very) degenerate (in the spirit of [Ivr8]). 

0.5 Ideas of the proofs and the structure of the paper 

To prove major theorems I follow scheme of [Ivrl, Ivr6]. Namely, first I replace operator A 
with its approximation A = A e (actually I consider two framing approximations) and then 
analyze the cut-off spectral trace 

(0.43) <f> e (hD t )xAt)Tr(e ih - ltA Q) 

with (p, x both supported in [—1, 1], % = 1 in [— |, |] and Xr(t) = x(t/T) etc; here Q is 
yU _1 /?,-pdo cut-off. 

Until 0.5.5 I will concentrate on the most difficult fi < eh~ 1 \ \ogh\~ 1 case. 
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0.5.1 Weak magnetic field approach 

First, I can introduce symbols = X{ — J2j Pij( x )d such that \{fi 2 a,qi}\ < C/i _1 where 
a = YlijS^PiPji Pi = Vi( x )- This means that evolve along Hamiltonian trajectories 
of u. 2 a + V with the speed not exceeding C/i -1 and for time T £ (1, Ti), T\ = efi, variation 
of qi and thus of Xi does not exceed Cu~ x T . 

Then I prove that variation of (some) linear combination of them, Yli^ilij * s exactly of 
magnitude /i _1 T. To do this I apply microhyperbolicity condition and it would be basically 
the end if I did not depend on r. Otherwise I need to consider an evolution of partial 
symbols /i 2 a^ (see (0.11)). The trouble is that |{/i 2 a, /x 2 a^}| is of magnitude 1 which means 
that only for time T < e in the classical evolution symbols /j 2 aj, are close to their original 
values and I can use microhyperbolicity. This would imply remainder estimate no better 
than 0(h l ~ d ). 

Remark 0.11. All above results for week, intermediate and strong magnetic field remain 
true if we replace partition 91 in Condition 0.1 by partition 971 and add Ch}~ d to right-hand 
expression of all estimates. 

To avoid this, I increase resonance groups to exclude intergroup third-order resonances; 
then |{/x 2 a, /x 2 a°}| are of magnitude /i -1 for (I, a) >z (2,0) and of magnitude log/i|~ CT 
otherwise. This means that /i 2 a° keep close to their original values for time T\ = efi as 
b (2,0) and T\ = efi l ~ 1 \ log/i|°" as (Z,er) -< (2,0). This explains where 01 came from 
and the first term in the right-hand expression in (0.26). 

Now, due to microhyperbolicity the shift in Yli ^'iVi ls °f magnitude exactly /i _1 T as 
T < 7\. In order to have it microlocally observable we need it to satisfy logarithmic 
uncertainty principle 

(0.44) ^T x e>Cu.- l h\\ogh\. 

Here e is a scale with respect to x (remember, that qi include £j via pj) and the right- 
hand expression is C/x — \ogh\ because we actually have /i _1 /i-pdos. So, for T £ [T ,Ti], 
T = Ce _1 /i| \ogh\ we can prove that 

(0.45) cj> e {hD t )xT{t)^ h ~ HA Q) 

is negligible for \ supported in [— 1, — |] U [|, 1]. Then construction of cut-off spectral trace 
(0.43) for T = Ti is reduced to its construction for T = T which requires analysis of the 
short-term evolution. 

On the other hand, due to rescaling arguments we know (0.43) modulo negligible for 
T st = e// -1 . If T st > Tq (which means exactly that e > Cfih\ log/i|), we are done. 
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Picking e = Cph\ \ogh\ we prove both estimates (0.26), (0.27) of theorem 0.3: the first 
term in the right-hand expressions is a remainder estimate Ch 1 ~ d T^ 1 for mollified operator 
and the second term is an approximation error estimate. 

This is done in section 1. 

0.5.2 Canonical form 

To consider a stronger magnetic field I need to reduce operator to a canonical form (kind of 
Birkhoff normal form). In contrast to 2-dimensional case resonances can appear and they 
are the real problem. Even in the smooth case removing "removable" fourth-order terms is 
not worth of trouble. Further, reduction involves Taylor decomposition with respect to p a 
(p — (pij • • ■ ,Pd)) and therefore each gained factor pj costs us differentiation with respect to 
x which in turn leads to division by e. Therefore in the heuristic sense I need \p\ < e. But 
in the microlocal sense \p\ x /i" 1 as fi < h~ 1 \ \ogh\~ l and therefore I must take e > /i" 1 
then. 

Note that in 2-dimensional case I was able to take smaller e = C(/i _1 /i| \ogh\)2 because 
instead of method of Taylor series I solved equation d^S = V — W along small circle (of 
radius x /i -1 ) where W is an average of V along this circle. In multidimensional case 
we have instead of circle a more general cyclotron orbit dcfij = fjdt which can be closed 
or everywhere dense depending on the ratio of fj] I was not able to repeat this trick in 
2r-dimensional case. 

Further, in the general non-smooth 2r-dimensional case without condition (0.22) I man- 
aged to make a reduction only in Ce- vicinity 4 ** and it is not enough even if due to weak 
magnetic field approach one needs only short-range propagation. 

Finally note that e = Cfi" 1 is better than e = Cph\ \ogh\ as fi > fi\. This explains 
that in order to derive sharp remainder estimate 0(/i~ 1 /i 1 ~ <f ) (for fx < /i" 1 ) I need in 
multidimensional case more regularity of coefficients than in 2-dimensional one. All the 
analysis beyond one of 0.5.1 should be done for \i > //* only. 

Actual construction is done in section 2 by /i -1 /i-Fourier integral operators in the smooth 
case and their "poor man's" equivalent in the mollified non-smooth case. 

0.5.3 Strong magnetic field: long range propagation 

I can use week field method to prove that for To < T < T\ with To = Ch\ \ogh\e~ l spectral 
trace (0.45) is negligible. In the intermediate magnetic field case it requires no modification. 

In the strong; magnetic field = eh 1 \ \ogh\ 1 < jjL < /ig = eh 1 one should take 

e = C(/j. — log/i|) 2 exactly as in 2-dimensional case. Then no modification is needed as 

4 ) Actually, one can increase it to e-vicinity as I > 1(d) but it is too high price to pay. 
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well with the only exception: for ffiR > 2 a very serious modification is needed in the part 
concerning localization of yU 2 a„. 

The problem is that y(/i 2 a° i , . . . ,/i 2 a° (n ) is not a proper symbol anymore (logarithmic 
uncertainty principle is violated) and this obstacle is not related to the smoothness. So, 
instead of "microlocal localization" (<£>(/i 2 a° i; . . . , /i 2 *^)) I consider "operator localiza- 
tion" y?(/i 2 a° i w , . . . , /x 2 a°™) and it takes some efforts to prove that for time t e [— 7\, Ti] this 
operator (applied to u) is negligible provided it was negligible at at t = 0; here and below 
b w is a Weyl quantization of symbol b. It requires to increase sightly e: 

(0.46) e > C^h] hgh\)h c ^ hlloghl . 

0.5.4 Short-range propagation 

Still, interval [—To, T ] is too wide. It would be nice to confine analysis to much smaller 
interval [-T *,T *] with T* = Ch\\ogh\ proving that (0.45) is negligible for T e [T *,T ] 
(actually one needs a bit more). Unfortunately it is not possible unless / is very large; I 
prove instead that (0.45) does not exceed Ch~ d (h/T) l \ log(h/T)\~ a . This statement implies 
immediately that for (I, a) y (1,1) 

(0.47) \F t ^ h -i T 4> e (hD t )xT(t) Tr(e ih ~ ltA Q) | < Ch l ~ d 

which immediately implies (due to Tauberian method) remainder estimate 0{h 1 ~ d T^ 1 ) for 
mollified operator in this case (but there will be still mollification error). 

0.5.5 Calculations 

So far, even if I recovered proper remainder estimates, the principal part of asymptotics is 
given by the same Tauberian method: 

(0.48) h- 1 J" (iWi T x T (t) Tr^^V)) dr 

with T = Ti or T = T due to negligibility of their difference. 

For weak magnetic field I had a standard expression due to standard results rescaled. 

In all other cases to get an explicit expression for the principal part of asymptotics I 
calculate u by a method of successive approximations on interval [— T ,T ] applied to the 
reduced operator. This is done in section 5. 
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0.5.6 Superstrong magnetic field 

This is completely different case. One needs to assume that = 1 here because arguments 
of 0.5.3 do not work anymore. Then the previous methods still work for [i < /^* = Ch' 1 . 

For /i > //*.* one needs to impose one of conditions (0.35), (0.36). So, I need to consider 
two different situations: 

(i) g 3 , Fjk are constant. Then making a reduction in e- vicinity and applying decompo- 
sition of the transformed u with respect to Hermite functions 

(0.49) U tlansiormed (x,y,t) = U aA X "' V" ifya (^h~^x') Vp {^ll'^y') 

ct£Z+ r 

we see that u a ^ satisfy a system 

(0.50) (hD t — B a (x", fi~ 1 hD")^u ajl 3 = b an {x" , ji^hD")^^ 

7 

with initial data 

(0.51) u a>p = 5 af3 5(x" - y") 

where B a = . 2(«j + ^f^fih — V w (actually one needs to transform fj and V first) and 
intertwining operators b a p are rather small. 

Similar system holds with respect to /3, y", t. 

Then unless 

(0.52) | Y^( 2a i + V fit** 1 + V\ < e 

j 

operators (hD t — B a (x", n,^ 1 hD"y\ are elliptic at energy levels r close to 0. Let 21 be a 
(finite) set of indices satisfying (0.52). Then expressing u a ^ with a ^ 21 via u a ^ with 
a e 21 and (doing the same for u a ^ with 21) I reduce the original system to a system 
for U = (u ai /3)a&&,pe2i which is analyzed in a standard way. 

(ii) Condition (0.36) holds. Then 21 = {a}. In this case I can reduce everything to just 
one equation. 

0.5.7 Some references and future plans 

The asymptotics with the principal part (0.4)* but without any remainder estimate was 
derived by G.Raikov in [Rail, Rai2, Rai3, Rai4, Rai5]. 
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Sharp spectral asymptotics in the superstrong magnetic field case (smooth) was derived 
by M.Dimassi in in [Dim]. 

In forthcoming papers I consider operators with magnetic field intensity matrices having 
incomplete constant ranks 2r. The typical example (already studied in [Ivrl, Ivr6]) is 3- 
dimensional magnetic Schrodinger operator. Canonic form of such operators contains "free" 
d — 2r-dimensional /^-differential Schrodinger operator. On one hand, this means that the 
remainder estimate could be only h times "better" than the principal part (while now it is 
yT x h times "better"). On the other hand, such operators require weaker smoothness and 
non- degeneracy assumptions. Under assumptions, similar to assumptions of this paper I 
am going to derive asymptotics with the principal part x {h + fi r h~ d+r ) and with the 
sharp remainder estimate 0(h l ~ d + [/h l ~ d+r ) as (I, a) y (1,2); further, under condition 
(0.35) for d — 2r > 3 and for d — 2r = 2, (/, a) y (2, 0) no non-degeneracy condition will be 
needed. Surely, only case d — 2r = 1 is generic. 

Later I hope to be able to consider (Fj^) of variable rank but remainder estimates will 
probably deteriorate. I also plan to consider Dirac operators with the strong magnetic field 
(as in Chapter 7 of [Ivrl]). 

1 Weak magnetic field 

As I mentioned, in the weak magnetic field we consider original (x, £) coordinates. As I 
explained we can take x-sc&le 



and therefore to satisfy logarithmic uncertainty principle we can take £-scale \x~ x . By no 
means this is larger than C{ji~ l h\ log /i|) 2 needed to consider /i _1 /i-pdo with symbol smooth 
in e-scale in both x and £ and therefore we cannot make a reduction here; only some kind 



(1.1) 



e > Cfih\\ogh\. 



of quasireduction 5 ) . 



1.1 Heuristics. I. Classical propagation 

Here I am interested in classical dynamics generated by symbol n 2 a° + V with 




5 )As in [Ivr6] in case the of weak magnetic field it is possible. 
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at energy levels close to for time t : |t| < T\ with T\ = e/i. On energy levels below cq we 
have 

(1.3) \ Vj \<Cn-\ 

At this heuristic stage I will not look at the smoothness (thus assuming that I is large 
enough) . 

First of al we can find j ]k G T l ' a such that for q 3 = f Xj — J2 k fjk{ x )Pk 

(1.4) {p m ,q J } = 0(fi- 2 ). 
Really, 

(1-5) {Pm,Pk} = H^Fkm 

and therefore {p m , q 3 } = ^ l 5 mj - /i" 1 J2j,k fjkFkm and one needs to take (f jk ) = {F jk )~ l . 
Then one can see easily that 

(1.6) {qj, q k } = {xj, q k } = fi^fjk mod 0(u~ 2 ) 
and therefore 

(1.7) /i 2 K, q 3 } = Oifi' 1 ), {V, q 3 } = Oifi' 1 ) 
which implies that 

(1.8) In classical evolution at energy levels below cq for time t : \t\ < efi increment of Qj 
does not exceed C/i" 1 ^! < 

(1.9) In classical evolution at energy levels below cq for time t : \t\ < e// increment of Xj 
does not exceed C/i~ 1 (|t| + 1) < e\. 

One can see easily that in the case of constant g^ k , Fj k 
(1.7)' /iVU} = 0, {V,q j } = u,- 1 Y,f j kd k V 

k 

and therefore 

(1.10) In the case of constant g jfc , Fj k the classical evolution for time t : \t\ < efi is approx- 
imately described by equation 6 ^ 

(i-ii) j t ^ = ^ l Y,f^ v 

k 

6 ^This equation describes approximately coordinates Xj but not velocities. 
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along which V is preserved. 

To describe evolution more precisely in the general case let us consider a point x. I 
remind that then VJt(x) is a partition of {1, ... , r} such that \(fi — fj)(x)\ < e iff i, j belong 
to the same element m of this partition. 

Then we can introduce 

(1.12) £?(z>0= a jk{ x )Pk(x,£,), j = l,...,r 

l<k<2r 

with complex coefficients ctjk G T l ' a such that 



(1.13) {QXk} = 0, {Cj,Cfc} = at S° d ^ f {(x,0,Pi = ■ ■ • = P2r = 0} 
where I remind that < means complex conjugation, 

(1.14) {C„ C k } = 0, {Cj, Cfc} = 2z/i- 1 5, fc mod 0(/i~ 2 ) as |p| < c^ 1 
and 

(1-15) «° = ^ E a iJfc (x)Cj(x,OCfc(^0 

with Hermitian matrices (ajk)j,kem close to a scalar matrices / m i# m of the same dimensions. 

Note that E° is a symplectic manifold and M d 3 x inherits its symplectic form, a volume 
form, and Hamiltonian field (we call it "Liouvillian field") respectively: 

(1.16) - ^Fjkdxj A dx k , 

(1.17) | detF|2 dxi A ■ ■ ■ A dxd = fx ■ ■ ■ f r ^fgdx, 

(1.18) £^ = -J2fAdrf)d k . 

Note also that 

(1.19) < d =^EM^c](x,oa(x,o 

j,k£m 

satisfy /x 2 {a°,a^} = 0(1), {F, aj,} = 0(1) so the increment of /i 2 a^ for time t does not 
exceed C\t\ and thus /i 2 a^ are moderate but not slow variables (and their increment for 
t — T\ could be rather large. 
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Now we can replace Q by 

(1-20) Q new = Cj ~ 5^ PjkmWVkVm 

kl 

with Pjkm £ T l ~ x ' a to make 

(1.21) {C„ Ck} = 0, {Cj, 0c} = 22/i" 1 ^-,. mod 0(/i~ 3 ) 

where for the sake of simplicity we use notations 772&-1 == Re f?2fc = ImCfc, & = 1, • • • , r - 
Then 

(1.22) /i 2 a° = -fx 2 ^ X M^KjCfc + g^ 2 J] a J« W 7 ?™ + • • • ' 

m _/',fcGm j,k,m 

(1.23) ^ = y(g) + ^a i (g)p i + ... 

where dots denote terms with 0(/i _1 ) gradients. 

These cubic terms in (1.22) and linear terms in (1.23) are the only sources of trouble 
because Poisson brackets /x 4 {a^, a^,} = 0(/i _1 ) and /x 2 {a^, V(q)} = 0(fi~ 1 ) where are 
defined by (1.19) with ajk(x) replaced by aj k (q). 

To get rid of these terms we can redefine Q again replacing them by Q + {S, Q} with 

(1-24) S = (1 X 'JkmnVkVmVn + ^ 7m??m- 

k,m,n m 

This replacement preserves (1.13), (1.14), (1.18) and modulo terms with 0(fi~ 1 ) gradients) 
it is equivalent to the replacement of fi 2 a + V by 

/i 2 a + V + {5, /i 2 a + V} = fi 2 a + V + /i 2 {S, a V}. 

We can always choose jj killing aj and if there is no third order resonances we can 
choose 7jfc m killing ctjkm] in the general case we can reduce /i 2 ao + V to 

(1.25) fi 2 J2a° m + n 2 ReJ2 X W?)CjCfcCn» + 

m n j,k,m£n 

where n denote third order resonance groups (see Introduction). 
Then, introducing 

(1-26) a°= f X^m+ Re X a ikm(q)(j(k(m 

mCn j,k,m£n 
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we see that /x 4 {a°, a°/} = 0(/i x ) and /x 2 {<2n, V(q)} = 0(/x x ) and therefore 

(1.27) In classical evolution at energy levels below cq for time t : \t\ < e/x increment of a n 
does not exceed C/i" 1 ^! < ei; 

(1.28) In classical evolution at energy levels below Co for time t : \t\ < e/x increment of 
a° = f Xlmcn a m does not exceed Cfi~ l (\t\ + 1) < ei. 

It follows from (1.6) that 

(1.6)' {q^q^^^fM+Oifi- 2 ) 

and then one can reintroduce qj d = Xj - J2 k fjk(%)Vk ~ Hk,mPikm{x)ri k ri m to make 

(1.29) {r ]k ,q j } = 0(v- 3 ). 
Then 

(1.30) {/xV + V, <&} = u 2 ^{a fem , + {V(q), g,} mod 0(/x~ 3 ) 

fc,m 

and microhyperbolicity condition (0.15) means exactly that 

(1.31) {/i 2 a° + V, Y^j PjQj} > ^oP' 1 as \x — x\ < e, |/x 2 a° — r n | < e and |/x 2 a° + V| < e. 
Thus, ^ ■ will really increase with the rate x /x" 1 for time \t\ < T\ — e/x. 

1.2 Heuristics. II. Classical propagation 

Now I would like to consider what happens for not very large smoothness. 
First of all, all our analysis obviously holds as (/, er) y (2, 0). 

Furthermore, as long as we do not need to analyze a° (i.e. if we have #91 = 1 in 
the microhyperbolicity condition) (^-smoothness is sufficient: to estimate by (^(/x" 1 ) the 
rate of change of q we need to have only {Cj,q k } = 0(/x~ 2 ) and {xj,q k } = (^(/x" 1 ). With 
microhyperbolicity we need {(j,q k } = o(/x~ 2 ) and calculation of {xj,q k } modulo o(/x _1 ). 
To do it we just need to mollify C-coefficients pj km with the mollification parameter o(l). 

If (/, a) z< (2, 0) and > 1 in the microhyperbolicity condition, we need to analyze 
evolution of a°. One can see easily that 

(1.32) {/x 2 a° + /, /x 2 a°} = 0(a l ' 1 \\ogfi\- (T ) 

because one can consider third-order and linear terms and replace them by their /x _1 - 
mollifications with respect to rj and Q; then for mollified terms this estimate would hold 
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and for mollification error it would hold as well since the gradient of the error would be 

log /ir). 

Therefore one can take 

(1.33) Ti = e^- 1 |log/i| <T 

and in the end of the day semiclassical error will be Ch 1 ~ d Ti 1 which is exactly the first 
term in estimate (0.26). 

1.3 Semiclassical propagation 

Now moving from classical dynamics to rigorous analysis I prove few statements, assuming 
that 

(1.34) Aio < < = log/i]" 1 , 

(1.35) e>Ch\\ogh\ 

where here and below /io, C are large enough constants, and e& are small enough constants. 
Proposition 1.1. 7 ) Let ji < eh~ l \ \ogh\~ 1 ^ and let 

(1.36) M> sup 2Vy% + e 

E jk 9 jk ^k+V=0 k 

with arbitrarily small constant e > 0. 

Let be supported in 5(0,1), 02 = 1 in -8(0,2), x be supported in [—1,1] 9 ^. Let 
T = Ch\logh\ <T<T[ = e . 

Then 

(1.37) \F t ^ h -i rX T{t){l - <p2 M T{x - x))<p 1M T{y - x)u{x,y,t)\ <Ch s Vr< ei 
where here and below t\ > is a small enough constant. 

Proposition 1.2. 10 ) Let condition (0.2) be fulfilled. Let 

(1.38) T = Ce- 1 h\ log h\ + C(fih\ log h\)^ <T<T 1 = e of i. 

7 ^Finite speed of propagation with respect to x; cf Proposition 2.1 [Ivr6]. 
8 ) Lower bound condition is not needed in this statement. 

9 'I remind that all such auxiliary functions are appropriate in the sense of [Brlvr]. 
10 ) Magnetic propagation; cf Proposition 2.2 (i) [Ivr6]. 
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Then 

(1.39) \F t ^ h -i T x T {t)(l - 02, m -ia/t(<?i - zi, ■ • -,Qd ~ Xd)Y x X 

u(x, y, t) (^ iM -imt(9i -x 1 ,...,q d - x d ))™ \ < Ch s W : \r\ < e x 

where here and below yU > 0, C\, M are large enough constants. Here and b w means Weyl 
quantization ofsymbolb and due to condition (1.38) logarithmic uncertainty principle 
holds and this quantization of the symbols involved is justified, B l means dual operator; 
(ii) In particular, for T > C\ this inequality holds for q x , . . . , q d replaced by x x , . . . ,x d . 

Proposition 1.3. n) Let v = #9T > 2. Let 



1.40) T x = ev-\ 



Cfi- 1 for (I, a) h (2,0), 

C/i 1 -'|log/i|- CT for (I, a) X (2, 0) 



and 

(1.41) T = Cv-\iih\\ogh\)^ < T x . 
Let <frj be functions as before. Then 

(1.42) \F t ^ h -i rX T{t)(l - 02,^(^0°! • ,/i 2 a°, - ^)) W x 

x u(x, y, t) (0i iJ/r (/j 2 a°i - n, . . . , /i 2 a°_ - 7v)) J*| < Vr, |r| < e x . 

All propositions 1.1-1.3 are proven by the same scheme as in [Ivr6]: 
To prove finite speed with respect to x one can use function 

(1-43) w (j^l + ^- Cq i 

with w function of the same type as used in theorem 3.1 of [Brlvr], q = ±1 depending on 
time direction and arbitrarily small constant e > 0. 

To prove that speed with respect to Q does not exceed C/i -1 one can use function 

(1-44) .((¥^'-4 

To prove that speed with respect to a„ does not exceed v given by (1.40) one can use 
function 

4 £iK -T n ,| 2 , 9A i * 



^ "( r +«■)»- c4) 

n )Evolution in , . . . , a°^; no 2-dimensional counter-part. 
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Proposition 1.4. Let (I, a) >z (1,2), microhyperbolicity condition be fulfilled and 

(1.46) n < n < €ih~*\ log/i| _ 2, 

(1.47) e = CV/i|log/i|, 

(1.48) T = e/i" 1 

where large enough constant C is chosen in the very last moment. Let us define 

(1.49) T\ = efj, if either (/, a) y (2, 0) or #91 = 1. 

(1.50) Ti = e/i'- 1 |log/i| CT i/(Z,<r) -< (2,0) and #9t > 1. 
T/ien /or T G [T , T x ], r G [-ei, ej 

(1.51) |F^ h -i T x T (t)rv«| < Ch s . 

Proof. (I) Let us first try function 

(1-52) w^(£,q-q) + e^) 

with I = £(y; 7~i, . . . , t„). One can check easily that it is an admissible symbol as long as 
T > T = e/i -1 , /i < 63/i _ 3| log/i| _ 3 and scales with respect to x, £ are 7 = e < ji~ l T and 
p = /i -2 < /i~ 1 T respectively; then fry > Cn~ l h\ \ogh\. 

Then our standard analysis (like in theorem 3.1 of [Brlvr], together with propositions 
1.1-1.3 imply estimate (1.51). 

(II) Now let us consider case e^h~^ | log h\~~. < fi < e\h~~- | log /i| _ 2 and e/i -1 < T < 
C(fih\ log/i|) 2 ; then 

(1.53) e > C(/x _1 /i|log/i|)3 > n~ x T, 

function (1.52) is not necessarily admissible symbol and one needs a bit more delicate 
arguments. 

Let us consider small vicinity of y. Without any loss of the generality one can assume 
that fjk{y) = as \j — k\ 7^ r (otherwise one can achieve it by a rotation). After this, 
without any loss of the generality one can assume that I = (£', 0) where x' = (x±, . . . , x r ), 
x" = (x r+1 , . . .,x 2r ) etc. 

12 )Weak magnetic field case. 
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I remind that due to definition qj = Xj — fjkPk + \ J2i k PjikPiPk, where (3jik are (mol- 
lified) coefficients chosen to eliminate linear with respect to (p±, . . . ,pd) terms in {pk,Qj}- 
Consider function 

(1.54) ™(^(l,q - q) + eq^j = w(^((£',x- x') -^ijfjkPk + ^^^jikPiPk)) 

j,k i,j,k 

and use a vector scale (7, p) with respect to (x, £) with 

(1.55) Pjlj — Cp~^h\ log/i|. 

Consider the first derivatives of the argument of to; we want them to be bounded after 
rescaling. We take 7, = 7, pj = p as j < r and 7, = p, pj = 7 as j > r with p > 7; then 

\fjk - fjk{y)\ < cp. 

Note that to have first derivatives with respect to x' bounded by Co7 _1 one needs to 
take 7 = p~ 1 T 13 ) as j < r. Further to have first derivatives with respect to £' bounded 
by Cop -1 one needs to take p = (p~ l T)z as j < r. Then first derivatives with respect to 
£" are also bounded by C^" 1 and to have first derivatives with respect to x" bounded by 
C Q p~ l one needs to have \djPk\ <cpasj<r,k = j + r; this could be achieved by gauge 
transformation. 

To satisfy (1.55) one needs then 

(1.56) T > Cp,{ir l h\ log/i|)3; 

one can see easily that T = T = ep~ x fits if and only if p satisfies (1.46). Then picking 
up 7 = mm(^Cph\ log h\, C/i _1 T^ and p = 7^ we see that (1.55) holds and (1.54) is an 
admissible symbol. 

Then our standard analysis (like in theorem 3.1 of [Brlvr], together with propositions 
1.1-1.3 imply estimate (1.51). □ 

Applying the same methods as above in the cases T > C{ph\ log/i|)3, T < C{ph\ log/i|)3 
respectively we can prove easily 

Proposition 1.5. 14 ^ Let (Z, cr) y (1,2), microhyperbolicity condition be fulfilled and 

(1.57) eth^l \ogh\~^ <p< eh~ l \ log/i]" 1 , 

(1.58) C'(p~ l h\ log/i|)3 < e < C'ph\ log/i|, 

(1.59) T = Ch\ \ogh\ ■ max(e _1 , (p^h\ log/i|)"^ 

13 ^Or less; we always take the maximal value. 
14 )Weak magnetic field approach. 
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where large enough constant C is chosen in the very last moment. Let us define T\ by 
(1.47) -(1.48). 

Then for T G [Tq,Tx] ; r G [—ex, ex] estimate {1-51) holds. 
Corollary 1.6. In frames of proposition l.J±, 1.5 

(1.60) |iWi T (x T (t) - XT (t))Tiju\ < Ch s 



.1 II 

2' 2J 



Vr : Irl < e 



where I remind that % is admissible function supported in [—1, 1] and equal 1 at 
1.4 Weak magnetic field: conclusion 

Since in frames of proposition 1.4 To = e^u" 1 and rescaling x 1— > /ix, £ 1— > t 1— > /it, 
/i 1 — > yu/i, T 1— > e leads to a standard Schrodiinger (i.e. with fi = 1), we can apply standard 
methods. Then 

(1.61) \F^ h -i T %r(t)rM < Chl ~ d ^ ■ M < e 

and applying Tauberian arguments we conclude immediately that 

(1.62) |r(^e)(r) - h' 1 [ (f^^xt^^u)) dr'\ < CT^h 

with any T G [To,Tx]. 

Now we need to calculate 

(1-63) h- 1 [ (F t ^ h -i T ,XT (t)T(iJu))dT'. 

J —00 

Using rescaling to standard Schrodinger operator we get that it is equal to 
(1.64) Yl Kn m h- d+2m+2n v 2n 

m,n>0 

(the rest of the coefficients vanishes). Taking only term n 00 h~ d we get standard Weyl 
expression J S w (T,x)ip(x) dx with an error 0(fi 2 h 2 ~ d ) which is 0([i~ l h l ~ d ) as \x < h~z. 
Taking only terms Koo^~ d + K io^ 2_ 'V 2 we commit error 0(h 2 ~ d +fi 4 h 4: ~ d ) which is 0([i~ 1 h l ~ d ) 
as n < h~5. One can calculate this correction term Kioh 2 ~ d fi 2 easily. Instead let us notice 
that k„ = / Xno(x,T)i(j(x) dx where Xio(x,t) depends only on g^ k ,V and dkVj calculated 
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at point x, and thus with indicated error to calculate (1.63) with T y instead of T and with 
ip = 1 and one can consider Schrodinger operator 

(1.65) A = J2g* k P j P k + V, gi k = gi k (y),V = V(y), 

Pj = hD 3 - VM - ^(«9^)(y)(* fc - y k ) 

k 

and for this operator expression in question is exactly £ MW (y, r). 
So, we have proven 

(1.66) |r(#)r) - J £ MW (x, r)i)(x) dx\ < CT{ 1 h 1 ~ d Vr : |r| < e 

where for ji < Ch~^ one can replace £ MW by £ w . Thus 

Corollary 1.7. Theorem 0.3 is proven under condition (I.46). 

The case of weak magnetic field is completed. In the case of intermediate magnetic field 
(as in proposition 1.5) I need to consider evolution until time T < eofih\ log/i| and according 
to proposition 1.2 evolution is confined to C/i _1 - vicinity of the original point. According to 
(weaken) proposition 1.3 evolution in (a^.)j = 1, . . . ,p is confined to ei-vicinity of r. Here 
we can use finer partition 971 rather than (because of the shorter interval). 

1.5 Propagation. Strong magnetic field 

Assume now that 

(1.67) eh^llogh]- 1 < /1 < C h~ l , 

(1.68) e>{Cii~ l h\\ogh\)^. 

In this case microlocally \p\ < C(fi~ 1 h\ log h\)^ rather than \p\ < C[i~ x even if e — 1; 
because of this statements of the previous subsection should be modified: 

Proposition 1.8. Under conditions (1. 67) -{1.6 '8) statement of proposition 1.1 remains 
true with M = C^ih\ \ogh\. 

Proposition 1.9. Under conditions (1.67)-(1.68) both statements (i),(ii) of proposition 
1.2 remain true. 
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So far T = Comuh\ log/i| and T\ = eji. 



If y u ^ = 1 the same arguments are sufficient to prove that XT{t)FuQ^ y = for 
Te [T ,Ti]: 

Proposition 1.10. Let = 1. T/ien under conditions (1.6f)-{1.68) statement of propo- 
sition 1.5 remains true. 

However, for v = #91 > 2 one needs to consider evolution of a n and the main trouble 
is that we can microlocalize fi 2 a n only in the box of size fi 2 e 2 x fih\ log/i| which is by no 
means smaller than e. We need instead operator localization: something like 

(1.69) Ilx,((/i 2 < - r ni , . . . , (/i 2 cC - T n „)^(x)u\\ < Ch s 

as \t\ < T where \ is our standard cut-off function and constant e could be taken sufficiently 
small (depending on the constant in microhyperbolicity condition). 

To do this we will need to consider something like x(B) with B a linear combination 
with positive coefficients of (/i 2 a™ — r ni , . . . , (/i 2 a™ — r n!/ ). Then 

(1.70) n~ 2 b(x, = f(x,p) = 2 S + 3 KkPiPjPk + + yu~ 2 6jPj-. 

Last two terms are barely in the principal part. 

We want to know how it evolves under propagator e~ l ^ h tA with \t\ < e 2 ■ Note that 
4>(B) = J c~ %tB (\){t) dz where is a Fourier transform and therefore we need to consider 
evolution of A under e~ trB propagation. 

This is much shorter evolution because instead of t : \t\ < e 2 we should consider r : 
M < C s \ log/i| since we can take |0(t)| < Ch s otherwise. Since B = yJn~ x ■ fih- fi~ 2 B where 
u~ 2 B is a standard /i _1 /i-pdo, we have a standard evolution now with just a bit larger time 
than the standard one: with time C\xh\ log/i| instead of C. 

Let us write down corresponding classical evolution for symbol h: 

df_ 

dp/ 

v- df f , df 



(1.71) 


d 


(1.72) 


d 


where 




(1.73) 




(1.74) 





opk ox 



J 
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and Ajk = Ylibik{Pk,Pj} is a matrix which is equivalent (in some metrics) to skew- 
symmetric and non-degenerate. 

Differentiating with respect to initial data z = (x,p)\ T= o we get 



(3<a v j 



/9/ v-. 



k k i 

where /c^ Q ^ = d^k and /5 < a means that [3j < ctj for all j and (3 ^ a. 
Condition 1.11. Let us assume that 15) 

(1.77) \d%b jk \ < C 1+|c V- |Q| | logep + C Va : \a\ < N = e | loge| 

and the same estimates hold for Fj k ; then this estimate with (/ — 1) instead of I holds for 

bijk- 

Further, let us assume that this estimate with I instead of I holds for b ; then this 
estimate with (I — 1) instead of I holds for b k . 

Under this condition one can prove easily that if 

(1.78) \p^\ + \x^\ < Cl +m e l ~W\ hge\- a + d 

holds for |/5 1 < \a\, then the second and third terms in the right-hand expressions of 
(1.75), (1.76) do not exceed 

(1.79) (7j +l/3| e ,+1 - |o| |loge|-' T + C r 2e 1 -*(b (a) | + |x (q) |) + C 2 

(actually 5 = here). 
15 ^This and the next assumptions can be arranged see [Brlvr] for the choice of mollification). 
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Remembering that |r| < e 5 and thus increment of x is less than e 1 s we see that if we 

replace Ajj by its value at some of the point of the trajectory, say by A^- = A^-| T=0 , then 
the error would not exceed (1.79) as well. Then since A generates a bounded group we see 
that inequality (1.78) holds for (3 = a too (because it holds for r = 0) 

So, by induction inequality (1.78) holds for |r| < e~ s and \/3\ < eo | log e | : 

(1.80) \dP<f) T \ < C\ +m e l - m \ \oge\~ a + d V|r| < e~ 5 V[/3| < e | loge| 
This implies immediately 

Lemma 1.12. (i) If estimate 

(1.81) \k ia) \ < M(c l 3 al e Hal \loge\- a + l) Va : |a| < e | loge| 

holds for symbol k then it holds for symbol k o <p T with |r| < e s as well. 

(ii) In particular, estimate (1.80) holds the principal symbol k® = k° o r of wannabe 

[r x h-pdo K = K(t) == e~ lTB K e lTB where K is fi~ 1 h-pdo with T l ' a symbol k°. 
Considering transport equations 

(1.82) ~d~T^ n ^ ^ l~"a-mkm = 

m<n 

we conclude that k — L\k§ satisfies 

(1.83) < MC 2 (4 Q '^" |Qhl |log£r ,7 + l) Va : |a| < e |loge| -2 
and solving (1.82) (with n — 1) and with zero initial conditions we get that 

(1.84) |A;i a) | < MC 2 \r\(c l 3 al e l -^- 2 \loge\-' T \r\ + l) Va : |a| < e |loge| -2. 

Then both symbols k = C\k\ and A; = £2^0 satisfy (1.84) with M replaced by C2M\t\e~ 2 
and then k 2 satisfies (1.84) with the same modification (and in both cases \a\ < eo| loge| —4) 
and so on: 

(1.85) \k^\ < MC 3 2n+|a| |r| n ^- 2n - |a| |log5|" CT + l) Ma : |a| < e | loge| - In. 
Then composing symbol 

(1.86) k ~ J^^hTK 
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we see that n-th term does not exceed (as 2n > I + 1) 

(1.87) M (^^)V|V +1 | loger 

and it becomes less than Ch s as n = n(e) == eo| loge| in the definition of e > Co(/i _1 /i| log h\)^ 
with large enough constant CV 
So in fact we need to take 

(1.88) k(r) d = f i^hTUr) 

n<n(e) 

and modulo negligible operator K(r) = k(r) w . Therefore we had proven 

Lemma 1.13. If symbol k° G T l,a and Kq = k 0w then for |r| < e~ 5 modulo negligible 
operator K = K{r) == e~ lrB K e tTB = k(r) w is [i^h-pdo with symbol k{r) defined by (l.< 

This lemma implies that 



;i.89) (j){B)k w = J ^r)k(fxhT) w e- i,lhTB dr 



where means Fourier transform of 0. 
Note that 



(1.90) [A, e ~ iflhTB ] = -ifih [ e-^ hT ' B [A, B]e~ iflh{T - T ' )B dr' 

Jo 

and that by construction [A, B] = /i _1 /i/c w with 

(1.91) k = W'V with k p G T l ~^ a . 

/3:|/3|<4 

Therefore 

(1.92) [A, e~ ifJlhl ~ B ] = —i/jh f k(/ihr') w dr' x e~ iflhTB 

Jo 

where k{r) is defined according to the above rules. 
One can check easily that 

(1.93) k(r) = M^'V 

a:\f3\<4 
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with 

(1.94) \k^(r)\ < C ,|Q|+1 |rr(^" 1 " |a| " 2n |log£r ,7 + l) Vn, a : |a| < e |loge| - In. 
Then one can prove easily that 

(1.95) \d™kf}{f)\ < M m C m+|Q|+1 ^- |a| - 1 |loger CT + l) 

Vm, a : 2m + \a\ < €q\ loge\ 

with M m = f sup J<m g J • j| m-5 ', == (\p\ + /x -1 ^ -1 <C 1 (due to our assumptions and choice 
of Co in e) and then one can prove easily that 

(1.96) M m < C m m\e' mlos ( losm+llosql ) . 
Similarly, one can prove easily that 

(1.97) |<9™fc2(r)| < M m C m+n+|a|+1 |rr(^" |Qh2n_1 | loge|- CT + l) 

Vm, n, /3 : 2m + n + < e | loge| 

(improvements are possible but pointless). 

Therefore, replacing fc njQ ,(r) it by Taylor decomposition with respect to r, one commits 
an error with an operator norm not exceeding 

(1.98) CM m (Cfxh) m h- mo ^- [ \T\ m \Ur)\dT 

m! J 

^ ^i Cm—m log (log m+| log — m\ log /z/i|+m log m 

where m is the number of the first dropped power and mo is some fixed exponent. To make 
the last expression less than Ch s it is enough to have 



m. 



(1.99) — mlog(logm + | loggQ — m| log/m| + m\ logm| < — C\ \ogh\ — C 

Replacing (logm + | logg|) by | logg| in the left side and minimizing it with respect to m 
we get m = m = Cf l {[ih)~ 1 \ \ogq\ and 

{yLh)~ 1 \\ogq\ > C\logh\ 

or 

I log ^ | > Cjjth\ \ogh\ 
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or 

(1.100) e > (\p\ + fi^)e c ^ loshl . 

Since \p\ < C(fi~ 1 h\ \ogh\)^ in microlocal sense we can rewrite (1.99) as 

(1.101) e > C(ii~ l h\ log/i|)^e C/ih|log/l1 

and here e varies from C/i _1 as p, = h~ x \ log/i| to h 1 ' 6 as /x/i = S\. 

In operator sense |p| < Cp~ x and we can rewrite (1.100) as e > Cp^e ^ 10 ^ but 
(modulo choice of a constant) it is exactly the same! 

1.6 Strong magnetic field. II 

As a result we get 

(1.102) <f,{B)A= ^{^h) m A m ^ m \B) 

0<m<rh 

with A m defined in a standard "smooth" way: 

(1.103) A m = Ad£ A; Ad° B A = A, Ad'£ A = [B, Ad^" 1 A] . 

The difference is that now 0( m )(.B) are still proper pdos but only with effective semiclassical 
parameter ph which is not very small. Anyway, symbols of A m satisfy 

(1.104) \a$\ < C m e-^M m < C m m\e-^\ logg|~ m 
Important thing is that 

Lemma 1.14. Let <p, <p\ be two fixed admissible functions with disjoint supports. Then 
(j)(B)A(j) 1 (B) = 0. 

Proof of this lemma is obvious. 

Now we can apply the same arguments as in Theorem 3.1 [Brlvr]. Consider zu^ 1 ^ — 
s) + Cp'H). Apply it to u ± = 6{t)u. Then 

(1.105) {hD t - A)u^ = Tih6(t)6{x - y) 
and applying /i _1 /i-pdo cutoff we get 

(1.106) (hD t - A)v ± = &hf(x, y), v = / = S(x - yW y (B)Q\ 
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where cp is supported in vicinity of A G M. while w is our standard function and then 



1.107) 



Reih' 1 



(hD t - A),w(e-\B - A) + C^H) 



v ± ,v ± ) \ < Ch s 



Transforming this inequality we get 



1.108) \{w'{e' 1 {B- A) +Cn~H) 



1A v ± ,v ± ) \ < 



Co 1 1( A,w (e" 1 (B - A) + Cofi-H) 



v , v ± 



Ch s 



where w' is derivative of w. 

Taking w to be a primitive of function ip 2 we get 

(1.109) ||¥>(e -1 (B - A) + C fi' 1 t)v ± \\ 2 < 

^( £ ^k m \\^(e-\B-\) + C^-H)v- 

^0<m<m 



Ch s 



with k = f C\fih\ logg)" 1 . 

Increasing C\ if needed (but neither Co depends on G\ nor C\ depends on Co) we get 
the same estimate with the right-hand expression 

(1.110) Cq" 1 (^k m \\^ m \e~ 1 (B-X) + C fi- 1 t)v ± \\) +Ch s . 

0<m<fh ^ 

Plugging ■^k m ip ( - m ^ instead of tp and (m — m) instead of m into modified inequality (1.108) 
we get the same estimate with the right-hand expression (1.109) for 

(—k m \\^ m \e-\B - A) + C /i- 1 t)^ ± | ' 
Vm! 



1.111^ 



Taking a sum with respect to m, < m < m, and again increasing C\ we get that expression 



(1.112) (— r m |^ (m) ( e_1 ( S - A ) +C /i" 1 t)w ± ||) 

0<m<m 

does not exceed itself multiplied by Cq C plus C/i s where C does not depend on Co- Taking 
Co large enough we get that expression (1.112) is negligible. So, 

(1.113) \H^\B - A) + C^~H)u^ y {B)Ql || < C/i s . 
Therefore we have proven 
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Lemma 1.15. Let ip be supported in e-vicinity of X and tp — 1 in 2e-vicinity of X. Then 
for T < Ti = Co-^min^,^- 1 ! \ogh\ a ) 

(1.114) \\(L-p(B))u$l(B)Ql\\<Ch s . 

Furthermore, one can plug B\ and B2 instead of B and A; then it follows from lemma 
1.15 that ip(B 1 )e lT2B2 (f(B 1 ) = as \r\ < e\ih~ x . Multiplying by y3 2 (r 2 ) dr 2 and integrating 
we get 

(1.115) <p(B 1 )ip 2 (B 2 )ip{B 1 ) = 
which allows us to apply simultaneous cut-off 

(1.116) $(£?!, . . . , B v ) A ^ <px(B x ) ■ <p 2 (B 2 ) ■ ■ ■ <p u (B u ) 
and we get 

Proposition 1.16. Let conditions (1.67), (1.101) be fulfilled. Let functions (pj be sup- 
ported in e-vicinities of Xj and let ipj = 1 in 2e-vicinities of Xj. Then for T < 7\ = 
emin(/i, log/i|°") 

(1.117) || (1 - $(£?!, . . . , B v ))u¥ y (B u B v )Ql\\ < Ch s . 

Consider now the proof that singularities leave diagonal - again as theorem 3.1 [Brlvr] 
and theorem 2.2.7 of [Ivrl] with = (£, x — y) etT~ l . We plug-in 

v = f u($(Bi, . . . , B U ))^Q^ localizing x near y and Bj near A.,- for all j, and take £ = £(y, t). I 
remind that Bj are linear combinations of a nk and Xj are corresponding linear combinations 

Of Tj. 

To recover the proof we need to estimate from below (^ 2 J2jkPj^jkPk w ^ w ) w hh w = 
zu((j)) w v, v = $(Bi, . . . , B v )u and Kjk = £(g^ k ); one can assume that V = —1. 

Due to proposition 1.16 one can replace w by $(-Bi, . . . ,B u )w with a bit more wide 
supports of ifj than those of <p~j. 

Let us introduce /i _1 -admissible partition ip u (x) and on each element of it apply gauge 
transformation making \Vj\ < C/i _1 . Rescaling x 1— »• fi(x — x v ) we get (/i/i)-pdos fipj and 
fi 2 J2j kPj-KjkPk while $(-B) is a legitimate (fih)-pdo as well. 

Then for fih < e\ 

$(£?!, . . . , B v f (/i 2 Pj K jkPk - eo + CfjLhj §(B U ...,B U ) 
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is non-negative operator and 

$(£?!, . . . , B„) V ^pJZ^p^S!, . . . , B v ) > l -e Q ${B u B V )*Q(B U . . . , B v ) 

in operator sense and therefore 

(1.118) ^ 2 J2pJ^w,w) > l eQ \\w\\ 2 -Ch s ; 

the rest of the proof needs no modifications. 

Now applying the same approach as in the proof of proposition 1.5 we arrive to 

Proposition 1.17. Let conditions (1.67), (1.101) and microhyperbolicity condition 0.1 with 
#91 > 2 be fulfilled. Then estimate (1-51) holds forT e [To,Ti] withT\ = emin(/!, log/i| CT ) 
and To = Ce~ l h\ \ogh\ < ei(fxh\ \ogh\)2 . 

2 Canonical form 

Now we need first reduce our operator to canonical form assuming that (/, a) y (2, 1) and 

(2.1) n\ = ei/T3|log/i|-3 < y, < p* 2 = eh~ l \\ogh\~ l 

(case jih > eh^\ \ogh\~ 1 and especially case jih > C we consider later). We assume that 

(2.2) e > Coji' 1 

rather than e = C(ji~ 1 h\ \ogh\) 2 as in [Ivr6]. This larger e makes certain parts of our 
construction much simpler but leads to a larger error. In this section we consider ji~ x h- 
quantization. 

2.1 Reducing "main" part 

First we need to reduce the main part. Consider point x in the vicinity of x and consider 
Hamiltonian fields H Vj of rjj reduced to E° where rjj are linear combinations of pj = £ 3 - — Vf 

V2k-i = f Re^fc, V2k = f I m Cfc; k = 1; • • • j r an d Ck are defined by (1.12)— (1.13) (rather than 
redefined by (1.20)-(1.21). 
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Consider symplectic map e tHq where q = q(x; 771, ... , r)d) is a quadratic form 

(2-3) q(x; ^, ...,%) = ^ <lij( x )ViVj 

hi 

with g y G J r/ ' ff . Note that 

(2.4) — (77* o e tH «) = {q, r] k } o e tH « = QijiVh Vk}Vi e tHq 

modulo quadratic form. Then using arguments of subsection 3.1 of [Ivr6] we can prove 
easily that 

(2-5) Vk °e H « = Y, Mx)rij + (3' kij {x, OViVj 

3 i,3 

with (3 jk G P~ u , P' kij G T [ - x!a . Here B = ((3 jk ) = e QA where A = (A jm ) with Q = (q kj ) and 
Aj m = 1 as j — 2% — 1, m — 2i, Aj m = —1 as j = 2i, m = 2% — 1 and Aj m = in all other 
cases. I remind that {r]j, 7?m}| s0 = ^jm- 

Obviously, for given matrix B G J rl,cr one can find symmetric matrix Q = (g^) of the 
same regularity such that B = e FA if and only if B is symplectic i.e. B^AB = A. In 
particular, one can transform (modulo quadratic form): rji into pi, 772 into P2P2 with /3 2 
disjoint from 16) , and 773 into p 3 + /3 3) iPi + ,#3,2^2, and 774 into (3 4 p 3 + Ai.iPi + /?4,2P2 etc: 

(2.6) B?7 m = T]' m = p m + ^ /SmjPj as m = 2k - 1, 

j<2k-2 

(2.7) Br] m = 7]' m = PmmPm. + ^ AnjPj as 77J = 2k 

j<2k-2 

where coefficients (3j k G T l ' a are chosen to make {r}j,r}' m } = Aj m . 

Now let us consider y u~ 1 /i-do Q = q w and transformation T(t) = e ltflh ("poor man's 
FIO"). 

Using arguments of subsection 3.1 of [Ivr6] one can prove easily that 

(2.8) w(-i) = ( V ' k r + (E^^) w 

i,j 

16 ) Since {pi,p k } is disjoint from for some k and we just rename x k into Xi and v. v. 
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modulo operator with the norm not exceeding Cfi 2 h. Then for 

(2.9) A = /i 2 (^6 jfc (x) W ) W 

jk 

we have 

(2.10) T(1)A°T(-1) = /i 2 (J2 b jk (x) V ' jV ' k ) W + /i 2 QT ^ ^) " 

modulo operator with the norm not exceeding Cfi~ l h where b'^ k G jF'~ 1,cr ; we will skip ' in 
what follows. 

Also, the same arguments show that 

(2.11) T(-l)V(x)T(l) = V(x) + YfflliY 

3 

modulo operator with norm not exceeding C[i~ l h + Cfi~ l \ logh\~ a where b" G T l ~ 1,c . 

I remind that while bj k , V are functions depending on x only, b'^ k , b" are complete 
symbols. In the end of the day we decompose them in the powers of 771 , . . . ,772,. as far as 
their smoothness allows. 

Due to construction 771 = Pi = £1 — Vi(x) 17 Note that we can get rid of V\ using 
gauge transform, after which condition (0.22) still holds. Exactly here we need condition 
(0.22) rather than more natural (0.22)'. Otherwise we would be forced to confine ourselves 
to Ce-vicinity of x which would cause problems in section 3. 

Then according to commutation property d xi r]2 = 1 on E°. Further, by construction 772 
does not depend on £1, £ 3 , . . . , £ d and therefore 

(2.12) rj 2 = a(x,^)(x! - A(a/,£ 2 )), «| E o = 1 

with a G JF Z,<T , A G JF' +1,<T and x' = (X2, ■ ■ ■ ,Xd). Let us redefine r] 2 = (x\ — A(a;', £2)). 

(2.13) Without any loss of the generality we assume that x = 0; otherwise in what follows 
one should replace x by x — x. 

Now we want to transform r\ 2 into X\. To achieve this goal let us consider T'(t) = 
e ifj,h {£i\) w _ Then we have the series of exact equalities 

(2.14) T / (-t)D 1 T(i) = D u T'(-t)X w T'(t) = A w , T' (—t)xiT' (t) =x x + t\ w . 
17 >I remind that r\k means T)' k denned by (2.6). 
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Therefore T"(l) will transform "new" 77 J into xi. Now we need to check how this transfor- 
mation will affect r)J with j > 3 and also a w , bj k , b*J k , and finally V and (b"r}j) w . 

Note that Hamiltonian function piA belongs to J rl+2 ' a 18 ). Therefore Hamiltonian map 
$t = e i/fp i A belongs to J rl+1 ' a as t G [0, 1] and then 

(2.15) \\T'(-t)r)jT(t) - ( Vj o d? f n < 

CC^" 1 ^) 2 ^ + e 7+2 " 3 | loge|"^) x (l +e 7+1 - 3 |loger") <C^ 2 h, 

(2.16) ||T'(-t)67 fc T(t)-(^o$ t ) w || < 

C(^ _1 /i) 2 (l +e l - 3 \\oge\-^ x (l +e I+1 - 3 |log£r") < C/i _1 /i, 

(2.17) \\T'(-t)bJ jk T(t)-(b ijk o^\\ < 

C(ii- l hf(l + e l - A \\oge\-~ a ^ x (l +e I+1 - 3 |log5r") < Ch, 

Furthermore, (2.16) holds for a instead of bjk- 

Therefore transformation of A given by (2.9) by T'(l) leads to the same expression 
(2.9) but all the symbols rjj, b'- k and b' i j k are replaced by r(- o $ 1; b'j k o $ x and o 
respectively and the total error does not exceed Cfi~ 1 h. 

Note that now 771, 772 are redefined as pi ° $1 = £ 1; p 2 o $ 1 = xi respectively. 

Let us consider rjj with j > 3. Note that {pi,f]j} = on E° and therefore 

(2.18) rjj o $ 4 = r)j + kjTji + ijikPiPk- 

k,i 

The same arguments hold for Xj with j > 2. Therefore, moving errors arising in the 
quadratic part of A into its cubic part, we can replace rjj o <3> t by rjj + fcj£>i (j > 3), Xj o $ t 
by Xj (j > 2) and in we can replace x\ o $! = xi — A(x', £2) just by — A(a;', £2)- 

From the beginning rjj with j > 3 could depend on ^ (see (2.6), (2.7)). After above 
transformation rjj o $ x could acquire "more" £1. Note however that originally {77^,772} = 
on E° and 772 o $ x = x\\ therefore {?7j o $ 1; xi} = on E° and in 77-,- o $ 1; which are linear 
combination of pt o $ 1; we can replace $1 by £& — \4(0, x'), k > 3 (modifying again 

18 )Really. A G JF'+ 1 - <T and multiplication by pi increases regularity by 1 due to inequalities 
\Pj \ < cjT 1 < e. 
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cubic terms); x' = (x 2 , ■ ■ ■ ,Xd)- Thus we arrive to the same quadratic expression but with 
Pi = £1, P2 = xi, p k = £ fe - 14(^2, • • • ,a?d;6), fcjfc = b jk (x 2 ,x 3 , . . .,x d ;£ 2 )- 

Repeating this process in the end of the day we arrive to the same quadratic expression 

with r] 2k -i = to = %k (k = 1, . . . , r), b jk = b jk (x", f") where now x" = (x r+1 , . . . , x 2r ), 



with 6 jfcm = b jkm (x,£) G J 7 ' 

Also, considering transformation of V(x) we see that at each step semiclassical error 
does not exceed C(fi~ 1 h) 2 s l ~ 3 \ \ogh\~ a < Cfi _1 h and in the end of the day we arrive to 



with W G P' a and bj = bj(x,£) G 
Therefore we proved 

Proposition 2.1. Let conditions (0.2), (2.1) and (2.2) be fulfilled, x be a fixed point and 
e < R < e 2 with small enough constant e 2 . Then there exists a bounded operator T such 
that 

(i) For operators Q — (f with symbol q supported in B(0, R) C M. 2d and Q' = q'™ with 
q[ supported in B[(x : 0) : 2C 1 R) C R 2d and equal 1 in B[(x : 0) : C x R) C M. 2d the following 
equalities hold modulo negligible operators 



£" — (Cr+l, • • • , £,2r)~- 



(2.19) 




(2.20) 




(2.21) 
(2.22) 



(I-Q[)TQ = 0, 
T*TQ = Q 



and modulo operators with norm not exceeding C\i 



(2.23) 
(2.24) 



T*i/jT ee V w , 
T*ATQ = AQ 



with 



(2.25) 




3 
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where bij G J rl,cr , bijk = bijk(x, £) G JF' 1,<J , 6 g JF /,ct , 6, = bi(x, £) G JF' 1,<J are real-valued, 
p2k-i = £k, P2k = %k (k = 1, • • • , r), Pi = pf, x' = Oi, ...,x r ),x" = (x r+1 , x 2r ) etc; I 
remind that ip = ip(x) is a smooth function; 

(ii) For operators Q' = q' w with symbol q' supported in -B((x, 0), Rj C M. 2d and Qi = qf 
with symbol q 1 supported in I?(0,2Ci-R) C R 2d and equal 1 in B(0,CiR) C M 2d 

(2.26) ||Q'T(1-Qi)|| <Ch s . 
Remark 2.2. (i) Calculations show that 

(2.27) 6y = ^o*o, b = VoV , 2jj = ijo^ 

where \& is a symplectic map M M D B(0,R) 0), 2Cii2) and ^ = 11**1 S o with 

S° = {x' = = 0} and 

(2.28) |detD*o| = /i---/r; 

(ii) I remind that one can rewrite quadratic part as 

(2.29) ^Y, Yl ^ lhD ") Z *j Z k -^hJ2 a kk {x", ^hD") 

me<M j,k£m k 

with Zfc = fi~ 1 hDf : — ixki Z k = H^hD^ + ^ x k] 

(iii) In particular, if #m = 1 for all m (no second-order resonances) one can rewrite 
quadratic form as 

(2.30) W'^^hD'^^Dl + ^xt). 

l<k<r 

2.2 Reducing next terms 

There is very little what we can do else because of all the resonances. We already finished 
transformation of the "main part" 

(2.31) Ao = /J 2 K( x "i ^hD'^PiPj + bo(x", ^ l hD") 

hi 

where one can rewrite quadratic part as (2.29); in what follows the quadratic part and the 
potential of operator generate corrections which could be written in the form 

(2.32) b" a {x",^ 1 hD"){^Y, 

a:3<H<r(7,ff)]+l 

(2.33) b , :(x"^- i hD")(p a r 

a:l<|a|<r(/,o-)l-l 
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modulo_ terms_ with norms not exceeding Cji l \\og ji\ a , Cji l \\ogji\ a respectively with 
b" a G b m e p-\a\,<r m Here and below \(l j(T )] = [/] un l eS s I G Z, a > in which 

case |~(7, a)] = 1 + 1 etc. 

The major problem are resonances. The second-order resonances (/; = fj as % ^ j) pre- 
vent us from diagonalization of quadratic form as in (2.30). We can only rewrite quadratic 
form as (2.29). 

Now it is easy to get rid of non-resonant terms of third-order in the main part and of 
first order terms in V. Let us consider transformation 



(2.34) U 



. il - 1 r 1 (^s 3 (/,ry,(')+Si(i",fV,e')) 



which is a pseudo-differential operator as /i > Ch s| log/i|s. As a result principal part Aq 
does not change, and the next terms 

(2.35) V=V<, «! = £ + e b '"(x",ap a 

o:|a|=3 o:|a|=l 

(which have norms 0(/x -1 )) are replaced by 

(2.36) A + K,/i 2 S 3 + Si} /W + ... 

with d° = Y2i jbij(x" ',$,") PiPj and "short" Poisson brackets {.,.}' (with respect to {x',C,') 

only). Here dots denote terms of the same types (2.32) with \a\ > 4 (modulo 0(fx~ l \ log /i| _cr )) 
and of type (2.33) with |a| > 2(modulo 0(ji~ l \ log/i|~°")). 

By an appropriate choice of S3 and Si we can eliminate all terms in the "main part" 
jj?a\ + {d°,fi 2 S3 + Si} save corresponding to third-order resonances; 

(2-37) Re( / u 2 EC jfc (^^"d0a) 

m,j,k 

where ( m = £ m + ix m , m = 1, . . . , r (we can assume without any loss of the generality that 
(2.38) {a , Cm}' = lfrn ^ n 

n 

where matrix (f mn ) has only positive eigenvalues if k . 

Note that sum (2.37) is restricted to m,j, k such that f m is not disjoint from fj + f k ; in 
particular b'C k = unless m, j, k G n G 91. 

One can get rid of non-resonant fourth-order (2.32) and second-order terms (2.33) in 
the same way but while the second-order resonant terms are only 
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(i) QCk with j, k belonging to the same m-group, 
the fourth-order resonant terms include 

(ii) CjCkCmCn and CjCKLCl witn fj not disjoint from f k + f m + fn and also 

(iii) CjCfcCmCn with fj + f + k not disjoint from f m + f n 

which we never intended to cover. The same is true for fifth-order terms too. 

Still we need to have some concern about only third, fourth and fifth order terms in 
(2.32) and second and third order terms in (2.33) because for fi > e/i" | log/i|~2 only these 
terms could be larger than Cfi~ 1 h. 

Thus we arrive to 

Proposition 2.3. Let conditions (0.2), (2.1) and (2.2) be fulfilled, x be a fixed point. Then 
there exists a bounded operator T such that 

(i) For operators Q = q w with symbol q supported in B(0, R) C M? d and Q' = q'™ with 
q[ supported in B((x, 0), l) C R 2d and equal 1 in B((x, 0), 1) C R 2d equalities (2.21), (2.22) 
hold modulo negligible operators and equalities (2.23), (2.24) hold modulo operators with 
norms not exceeding C(/i _1 /i + \ log/i| -fJ ) with 



with A defined by (2.31) (and equal to (2.29)) and with symbols 0^3,0'^ E P |a| l/3|,,J and 



(2.39) A = A + n 2 Re 




b'^(x",^hD") (CCT+ 



Re 





unless 



(2.40) 



with arbitrarily 
(ii) Further, 



small constant e > 0. 



(2.41) 





a,/?:l<|a|+|/9|<r(«,<7)l-l 




a 




Tl<H<|a|+|/9| 

(iii) Statement (ii) of proposition 2.1 remains true. 
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I remind that only cubic terms are obstacles in the proof that does not change much; 
we prove this statement only for a°. 

2.3 Strong magnetic field case 

The same construction works in the strong magnetic field case 

(2.44) fi* 2 = eh-^logh^ 1 <fi< iil = Ch~ l 
but in this case one needs to take 

(2.45) e = C(/i- 1 /i|log^|)3 
(or larger) and restriction to Pj is now 

(2.46) \pj\ < Ce. 
However, in operator rather than microlocal sense we have still 

(2.47) IbjH < C/i- 1 

on energy levels below cq and we should estimate perturbations based on this estimate 
rather than (2.46). 

Note, that now we need to take in account only the third and fourth order terms in the 
first sum in the right-hand expression of (2.39) and only the first and second order terms in 
the second sum in the right-hand expression of (2.39). Similarly in (2.41) we need to sum 
as \a\ + \/3\ < 2. 

2.4 Superstrong magnetic field case 

The same construction also works in the superstrong magnetic field case 

(2.48) fi>f4 = eh~ l 

we should take again e by (2.45) and inequality (2.46) again holds in the microlocal sense 
but (2.47) is replaced by 

(2.49) Ib7|| <? = cvx-3/ii 
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Really, due to our assumption (0.34) 

(2.50) A = fi 2 M*"> ^ hD ") Z j Z k + b' (x", n^hD") 

mem j,kem 

with symbol b' e JF' ,<T . 

Then the principal part of operator A is of magnitude Cfih while norms of perturbations 
Ai, A 2 , A 3 , A 4 do not exceed Cfihq, Cfih<? = Cti 2 , C/i/i? 3 , Cfih^ 4 = C/i" 1 ^ 3 respectively 
and only the last one could be ignored for really large fi. This is the reason beyond correction 
in theorem 0.9. It is an also one of the reasons beyond last term in the right-hand 
expression of estimate (0.39) and the corresponding extra-smoothness requirement there to 
get the best possible estimate; another reason is mollification error. 

Anyway, in both strong and superstrong cases we arrive to 

Proposition 2.4. Statements of proposition 2.3 remains true in both strong and superstrong 
magnetic field cases with (2.24) holding modulo operator 1Z such that 

(2.51) \\Kv\\ < c(e l \ log/i|- CT + fihe l \ log/il"' 7 ) ||v|| + e~ l \ \ogfi\~' 7 \\Av\\ \fv 
with e given by (2.45). 

2.5 Canonical form. Special case 

One can consider it as a pilot-model for this section. 

If gi k , Fjk are constant then Vj(x) are linear functions and A is transformed into exactly 

(2.52) J2 fi(h 2 D> + ^) 

l<j<r 

by /i _1 /i-metaplectic transformation which consists of the following steps: 

(i) Change of variables (x,fi~ x hD) i— > (Qx,Q t ~ 1 fi~ 1 hD) transforming gi k into Sjk and 
F into matrix with Fjj +r = fj, Fj +r j = —fj and other elements 0. It transforms V(x) into 
V(Qx). 

(ii) Gauge transformation (multiplication by e lflh s ^ with quadratic form S(x), trans- 
forming Aj(x) into and Aj +r (x) into fjXj for j = 1, . . . ,r. Then A is transformed into 

(2.53) h 2 (d 2 + (D j+r - fjUh^Xif) 

l<j<r 
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and V(x) is preserved. 

(iii) Fourier transform: (x', x"\ fi~ 1 hD', fi~ 1 hD') i— > (a/, —fi~ 1 hD'; n~ 1 hD', x) transform- 
ing A into 

(2.54) h 2 ( D J + ^ h ~ 2 i x J+r ~ fi*i) 2 ) 

and 7(x) into V(x', -fi^hD'). 

(iv) Change of variables (x', x") /i^hD', ji^hD') i— > (x' — Kx" ,x'\ ji^hD' , ji^hD" — 
K t fi~ l hD') with Kjj +r = fj and other elements transforming A into 

(2.55) £ (/^J + /| A?) 

and 1/ into V(x' — Kx" , fi~ 1 hD" — K t fi~ 1 hD') with Weyl quantization. All other steps are 
the same as above. Moreover, if we do not look for uniformity, we can replace e by in 
condition (2.40). 

3 Mid-range propagation 

By mid-range propagation we assume propagation for A and for T*AT as T G [Tq, T[] with 

(3.1) T{ = C( f ih\\ogh\)K 

(3.2) Tq = Ce _1 /i| log/i| < Cfj,h\ logh\; 

we are interested almost exclusively in extending proposition 1.4. Here T[ is what used to 
be To and Tq is a natural improvement of To as if A was // -1 /i-pdo A(x", fi~ 1 hD") satisfying 
microhyperbolicity condition. 

However Tq is still pretty large unless we increase e and thus smoothness assumptions. 
Instead in the next section we consider (smoothness-dependent) rate of decay of the left- 
hand expression of (1.51) on interval [Tq^T"] with T" = Tq and Tq" = Ch. 

3.1 Intermediate magnetic field case 

Let us consider intermediate magnetic field case (2.1). 

Consider a point (#,£)• During time T\ point (x",£") will stay in C([i~ l h\ log h\) 2- 
vicinity of (x , £ ) in both classical and microlocal senses (see proposition 1.2 for the proof 
in the microlocal sense). 
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According to subsection 2.2 

(3.3) A T = T*AT = A + TZ, 
with 

(3.4) A= Am ' 

0<m<l(l,a)~\-l 

(3.5) 11 = p{x,iT l hD), p E fj,- l \\ogfx\- a J r0fi , 

Ao given by (2.31) (where one can rewrite quadratic part as (2.29) in the general case and 
as (2.30) as #m = 1 for all m) and 

(3.6) A m = p 2 Re b'^(x",^ 1 hD")(CC W r+ 

Q,/3:3<|a| + |/3|=m+2 

a,P:l<\a\+\P\=m 

(see (2.39)) and Q = ^ + iXj, Q = £j - ixj, j = 1, . . . , r. 

Remark 3.1. Surely we need to remember that A and Aq- = T*AT are close but not 
equal. In [Ivr6] we considered A and took TAT* as approximation 19 ) but now we have 
the global reduction only in the special case of condition (0.22) while until next section I 
need only (0.22)'. The following statements hold both for U' = XJr = T*e~ h At T and its 
Schwartz kernel vl = uq- and for U' = e lh and its Schwartz kernel vl . 

Proposition 3.2. Let u' be defined in remark 3.1. Let <f>\ be supported in 5(0, 1), 02 = 1 
in 5(0,2), x be supported in [—1,1]. Let Qk = <j>k,Mp- 1 T(%" — % , £,x — £ ) w - Then for 
To == C(fih\ log/z)5 < T < T[ = CqTq and large enough constant M 

(3.7) \F t ^ h -i TXT (t) (1 - Q2x)u aP {x, y, t)Q{J < Cfi~ s Vr < c. 

The same estimate holds if we replace 4>*(x") by cj)*(x'" , fi~ 1 hDiv) with x" = (x"';x IV ). 

Proposition 3.3. Let u' be defined in remark 3.1. Let m = #9Jt > 2 and let <pj be functions 
as before. Then for C(fih \ log/i|)a <T<€q 

(3.8) \F t ^ h -i rX T(t)(l - 2 ,r(/i 2 < - ri, . . ,,/i 2 al m - r m ))"x 

x u'(x, y, t) (<f) hT (n 2 a° mi - r 1; . . . , /i 2 a^ - r m ))™'| < C^~ s Vr, |r| < e 1 . 

19 >1 remind that actually we consider A e instead of A, so we already have an approximation. 
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Both propositions 3.2-3.3 are proven by the same scheme as in [Ivr6]: 
To prove proposition 3.2 one can use function 

(3-9) T2 ' +^-^| 

with function of the same type as used in theorem 3.1 of [Brlvr], q — ±1 depending on 
time direction and arbitrarily small constant e > 0. 
To prove proposition 3.2 one can use function 

( 3,o, ^( ^';r Tm/ +^- c 4 

Because T < eo and we picked up v = 1 rather than given by (1.40), third-order resonances 
do not pose any problem here. Even if at this stage decomposition in Hermitian functions 
is not very useful, we can study propagation using only functions of (x", t). 

Proposition 3.4. Let u' be defined in remark 3.1. Let microhyperbolicity condition 0.1 be 
fulfilled. Then for T G [Tq, T[], t G [-e', e'] 

(3.11) \F t ^ h -, rXT {t)T"Qu\ < Cfi~ s 

where Q = Q(x, fi~ 1 hD) and here and below V" is a partial trace (with respect to x" only). 

We can assume that £(x , |"; f±, . . . , f m ) has zero ^''-component and is equal to (£', 0). 
Otherwise one can achieve it by a linear symplectic transformation in (x",£") and corre- 
sponding FIO (which will be a metaplectic operator in this case). Then using 

(3.12) wLT-\i',x" - y") ± etT- 1 ^ 

which is admissible pdo-symbol one can prove proposition 3.4 in the standard way. Since 
proposition 3.3 holds with with DJl instead of 9T, the same is true for 3.4 as well (but not 
for corollary below). 

Propositions 1.5 and 3.4 imply immediately 

Corollary 3.5. (i) Under microhyperbolicity condition 0.1 estimate (1.51) holds for T G 
[Tq,Ti] with T\ defined in section 1: T± = efi as = 1 and T = e G (/x, logfi\ a ) as 
#9T>1. 

(ii) Also estimate (1.60) holds with T replaced by Tq. 
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Corollary 3.6. Theorem 0.1 holds in the case of intermediate magnetic field (and thus is 
proven). 

Proof. Picking e = Cfih\ \ogh\ we get approximation error Ch~ d (/j,h) l \ logh\ l ~ a and Tq = 
eyU" 1 . Then we again have estimates (1.61), (1.62) and need just calculate expression 
(1.63) which is equivalent (modulo negligible) to (1.64) where without spoiling remainder 
estimate we can take only terms with m = and n < n(5) as fi < h~ 1+s with arbitrarily 
small exponent 5 > 0. The rest repeats arguments of subsection 1.4. □ 

3.2 Strong magnetic field case 

Let now eh~ l \ \ogh\~ 1 < fi < e\h~ l . 

Then analysis makes real sense only for e ^> ([i~ l h\ log h\) 2 because T[ x Tq for 
e = (/i _1 /i| log/i|) 2 . Still even in the latter case we need to state propositions similar 
to propositions 3.3-3.4. 

Proposition 3.7. Proposition 3.2 remains true with Tq,T[ defined by (3.1), (3.2) as 



with small enough constant e\ > 0. 

Proposition 3.9. Proposition 3.4 remains true with Tq,T[ defined by (3.1), (3.2) under 
assumption (3.13) as #91 = 1 in microhyperbolicity condition 0.1 and under assumption 
(3.14) as #01 > 1 in microhyperbolicity condition 0.1. 

Proofs of propositions 3.7, 3.9 repeat those of propositions 3.2, 3.4. Proof of proposition 
3.8 repeats arguments of subsections 1.5, 1.6 dedicated to the proof of proposition 1.17. 

3.3 Superstrong magnetic field case 

This case is split into two subcases: 



In the former subcase we have proposition 3.7 and its long-range version proposition 1.9 but 
we neither have nor can have (in general) proposition 3.8; further, proposition 3.9 is proven 
(as #01 = 1) under too restrictive condition 0.1 rather than more appropriate condition 
0.2. In the latter subcase we have nothing about propagation at all. 




C h~ l . 

- T[ defined by (3.1) as 
e 1 h~ 1 



(3.15) 
(3.16) 



/ig = e\h 1 </i < n* A = Coh 



-1 



-1 
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3.3.1 Hermitian decomposition 

In both cases, however, we can apply very usefully decomposition (0.49) 
(3.17) ur(x,y,t)= ^ u aP (x",y",t)T a (x')T p (y') 



with T a (x') = v a (fi2h ax') and T a (x) = Y[i<j< r v aj( x j) an d Hermite functions v aj . 
Note that 

(3.18) A (v a (x")T a (x') 



/3:\/3\=\a\,\/3-a\=0,2 

where 

^^(2«j + l)bjj as (3 = a, 



(3.19) A 



o 

0a 



2y / a j f3 k b jk as $ - at = 5 ik - 

Further, H n =' {^2 a .\ a \ =n v a (x") r Y a (x')} are invariant subspaces of Ao and 

(3.20) Spec(A° af3 (x",e)) a « = + H = 4 

lnl n k _ J 

3.3.2 Special case of constant gi k , Fj k 

This is far the easiest case. First, we can assume with no loss of the generality that 

(3.21) A = A . 

Really, no cubic terms appear from "kinetic" part of Hamiltonian ^ . g^ k PjP k and first- 
order terms appearing from potential are eliminated in the process of reduction. Further, 
we can skip an error 0(e'|log£|- CT ) and for (I, a) y (2,0) we can skip quadratic terms 
appearing from potential because in operator sense they are less than fi~ 2 (A° + 1) and 
since our reduction in this case is global, we can just take approximate operator as TAqT*. 
Further, we can assume with no loss of the generality that bj k = fj5j k an d then 

(3.22) A{v a (x")T a {x')^ = (J2( 2a i + l )fi^ + qo)v a {x") ■ T a {x'). 

j 

Now we are dealing with r-dimensional scalar y u~ 1 /?,-pdos and the proofs of the following 
two statements are obvious: 
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Proposition 3.10. Let g^ k , Fjk be constant and one of conditions (3.15), (3.16) be fulfilled. 
Let Q = q(x" , ii~ x hD"). Assume that on support q ellipticity condition (0.41) is fulfilled. 
Then for T > Ch\ log h\ , r £ [— e', e'] (with small enough constant e' > 0) 

(3.23) \F t ^ h -i T {xT(-t)urQi)\ < CTpT 5 . 

The proof is based on standard elliptic arguments. 

Proposition 3.11. Let g^ k , Fjk be constant and one of conditions (3.15), (3.16) be fulfilled. 
Let Q = q(x" , fi~ 1 hD") . Assume that on support q microhyperbolicity condition (0.19) is 
fulfilled. Then for T £ [Tq, T x ], t £ [-e', e'] 

(3.24) \F^ h -i T (xtWutQI) I < Cfi~ s . 

The proof is based on our standard microhyperbolicity arguments. 

3.3.3 Not-ultrastrong magnetic field case 

Assume that of condition (3.15) holds. 

Consider operator A? as yU~ 1 /i-pdo with respect to x' with "matrix" symbol with values 
in £(H, H) where H = L 2 (W X ,) and £(Hi, H 2 ) is the space of linear operators from M 2 to Hi. 
Obviously A falls into this class, but we can also replace At by Aq + Q*(At — Ao)Qi where 
Qi = I — (^(e^x', e~ 1 fi~ 1 hD f ) and is supported in 2Ci-vicinity and equal 1 in Ci-vicinity 
of £ W^,,^,,. Applied to u, Qi produces negligible output (after applying F t _> h -i r XT(t) 
with r < c and Ch\ \ogh\ < T < chr s . 

Then standard elliptic arguments imply 

Proposition 3.12. Let condition (3.15) be fulfilled. Let Q = q(x" , fi~ l hD") . Assume that 
on support q ellipticity condition (O.4I) is fulfilled. Then for T > Ch\ \ogh\, r £ [— e',e'] 
(with small enough constant e' > 0) estimate (3.23) holds. 

Furthermore 

Proposition 3.13. Let condition (3.15) be fulfilled. Let Q = q(x" , fi~ 1 hD") . Assume that 
on support q microhyperbolicity condition 0.2 is fulfilled. Then for T £ [Tg, 7\] , r £ [— e', e'] 
estimate (3.24) holds. 

Proof is based on standard arguments of [Brlvr] combined with observation that the 
required microhyperbolicity condition 

(3.25) ((£A (x",£'))v,v) > e \v\ 2 - C\(A - r)v\ 2 W £ M 
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is equivalent to condition 0.2 (with £ depending on (x", £") only). Here (•, •) and | • | denote 
inner product and norm in HI. 

Note that we need the single direction t. We just cannot localize in (/i 2 a° l5 . . . , /U 2 ^ 
even after reduction. This is really frustrating because Poisson brackets { fJ?a° n . , j^ 2 a° nk } are 
0([i 3 h\p\ 4 ) which is 0(fih 3 \ \ogh\ 2 ) even in microlocal sense. 



3.3.4 Ultrastrong magnetic field case 



Assume that of condition (3.15) holds. Then we should take into account that unless Fj^, g 
are constant (the case we already considered), generally variations of fJ,hfj are of magnitude 
fih. To overcome this difficulty we need to assume that condition (0.36)-(0.37) holds. 

We still need to deal with the fact that we have only local reduction and thus Ar differs 
from A by operator A! with symbol belonging to cuJF ' with 

(3.26) u = e l \ \ogfi\- a + /ihe l \ log/x|" CT . 

Lemma 3.14. Let b G J 70,0 be supported in {\x'\ + |£'| < ce}. Then operator 

(3.27) blp : K 3 v -> {b w (v <g> Tp), T a ) G K, K d = L 2 (Rl„) 
has norm 

(3.28) ||6*j|| < Ce" e|Q - /31 . 

Proof. Obviously, it is sufficient to prove lemma as x' G M 1 (with x" G W). Then since 

(3.29) Hf£ p = (2(3j + l)^KT p , H 3 = /jLh(\Q\ 2 ) w = \?x) + h 2 D 2 
we conclude that 

(3.30) ^K-^rMioi 2 ,^ 

and continuing this process we conclude that b™p = (atj — Pj)~ n (b( n ))™p with 6( ) = b and 
6( n ) = {|Cj| 2 , &( n -i)}- One can see easily b( n ) G C n n!jF ' uniformly with respect to n 
20 > operator norm of bjs in £(H <g> K, EI <g> K) does not exceed C n n! and then operator 
norm of 6™^ in £(K, K) does not exceed \aj — (3j\~ n C n n\. It reaches minimum e~ e \ a i~$i\ as 
n = C~ l \aj - (3j\. □ 

20 )With simultaneous decay of the number of derivatives checked, but everything works as long as n < 
e|loge|. 
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So far we did not use condition (0.36)-(0.37). 
(i) Assume first that all fj are disjoint. Then 

(3.31) Ao= b j (x", / 2~ 1 hD")(fi 2 x 2 J +h 2 D^ + b (x",i2~ 1 hD") 
and 

(3.32) Ao(v(x")T a (x')) = (W a v)T a (x f ) 
with 

(3.33) W a = (V + ^2(2aj + I) ft) o tf . 

Proposition 3.15. Let conditions (3.16) and (0.36) — (0.37) be fulfilled (with fixed a). Let 
all fj be disjoint. Let T = fi m with large enough m. Then 

(i) As (a, (3) 7^ (a, a) 

(3.34) F t ^ h -i T x T (t)u a p = E ax F t ^ h -i T x T (t)uaaE t f3y 

Vr : |r| < e ah in Q x Q 

where Q = B(z ,Cqs) C is a domain in which reduction is done, z = (x ,x ), 

E a = E a (r, x" , ji^hD") is pdo with norm 

(3.35) \\E a \\ < C(^/0~ 1 ((/4~ 1 # |a ~ a| + e~ |a_s| w). 

(ii) Further, u m satisfies 

(3.36) F t ^ h -i T x T (t)(hD t - W 5t (x\^ 1 hD") - W a (x" \ ^hD'^u** = 

Vr : |t| < e'/i/i in f2 x f2 

and adjoint equation with respect to y" where W a is defined by (3.33) and W' a is pdo with 
symbol belonging to (uj + /i 2 )jF°'°. 

Proof. Plugging representation (3.17) into equation (hD t — A-r^u = and using (3.32)- 
(3.33), form of A and applying lemma 3.14 to Aq — A we arrive to system 

(3.37) (hD t -W a - W a ) u a * = B ^ « 7 * in ^ x R 2r 
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with W' a not including hD t . Here and below all operators are acting with respect to x" . 
We also get an adjoint system 

(3.38) u*/3 (hD t - W/3 - W'X = u *i B iP in R2r x Q 

with operators acting with respect to y" . 

Applying F t ^, h -i T XT(t) with T described above, we arrive to system 

(3.39) (r-W a - W^v a * = ^ B ai in ft x M 2r , v a(3 = F t ^ h -i T x T {t)v aP 

and similar adjoint system. 

Due to condition (0.36)-(0.37) for |r| < e[ih all operators (t — W a — with a ^ a 

are elliptic with inverse operator norms of magnitude (|a — alfih)^ 1 . 

Note that even in the worst case B a p are operators with symbols belonging to [ih([i~ l h\a\) 2 JF ( 
Then by successive approximations v a * = T a Va* asa^a where E a is the sum of terms of 
the following type: 

(3.40) (r - W a y l B' aai {r - W^B'^ ■ ■ • (r - W^B'^ 

with k > where -B^ are operators with symbols belonging to ^h(^" l h\a\)^ J^ 0,1 . Similarly, 
v*p = v*a Tt. For the sake of simplicity we included W' a into W a . Statement (i) follows 
then from these two formulae. 

Plugging (3.33), (3.17) into equation {hD t —Aq- x )u = and adjoint equation with respect 
to y, we arrive to equation similar to (3.36) (as well as adjoint equation with respect to y): 

(3.41) (t- Wa)va* = WaVa* m H X f 2r 

where Wa is a sum of terms of of (3.40) type with an extra factor on their left: 

(3.42) (^h^BUr ~ W^B'^r - W^B'^ ■ • • (r - W^B^ 
with B' aa of the same type as above; p = here. Note that 

(3.43) (r- W a )~ l = (W a - W^- 1 - (r -Wa)- 1 ^ -W^iWa-Wa)" 1 in ft x R 2r . 

We apply this formula to (r — W aj J -1 . Then we drag (r — W a ) to the right. If it dies 
at commuting, we get the same expression with the same k but with p replaced by p + 1 
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because in commuting we gain e 1 fi h factor. Continuing this process with (r — W ak _ 1 ) 
etc we arrive to 

(3.44) (t - Wa) v a * = (W'l + (r - W s ) + W^j v a * in Q x M 2r 

where W^' does not depend on r, W^" is the sum of products of (3.42) type with p replaced 
by p + 1, is the sum of products of (3.42) type with p replaced by p + 2 and k replaced 
by k — 1. Multiplying by (I + W&) one can rewrite this equation as (3.41) with W a replaced 
by W a + W a ' and W^' of the same type with p replaced by p + 1. 

Continuing this process we arrive to negligible W'£. Statement (ii) is proven. □ 

Remark 3.16. It follows from construction that modulo (u + n~ l h)T®'° added correction 
is equal to 

(3.45) W' a = BUWa-W^B^ 

\a—a\=l 

with B a p appearing exclusively from cubic terms in A. Further, note that these terms are 
b m ;jk Z m Z j Z k an d Ki-,jk Z m Z *j z l with fm not disjoint from fj + f k . Condition (0.36)-(0.37) 
implies that then a m = and either otj = or a k = 0; therefore cubic terms, applied to 
v(x")T a (x'), produce 0, and therefore W a G (cu + /i -1 /i).F ' . 

However, correction W a = h 2 Y^j K j Z j Z j m °d 0(u + fi^h) is already generated in the 
process of reduction in section 2. Still we see that it in the case of constant g^ k , Fj k . 

(ii) Consider now a more general case when some of fj are not disjoint; assume that fj 
is disjoint from the rest of eigenvalues as j = 1, . . . , p and is not as j = p + 1, . . . , r. Then 

(3.46) A = ^bj{x"^~ 1 hD"){^x 2 j + h 2 D]) + 

i<j<P 

b jk (x", fi^hD") (Z*Z k + fih5 jk ) + b (x", fi^hD") 

p+l<j,k<r 

Note that condition (0.36)-(0.37) yields that otj = as j > p+ 1 and thus Aq restricted 
to v(x")T a (x') equals 

(3.47) A = bj(x",/i- 1 hD")(2a j + l); +b {x'\ fi^hD") 

i<j<P 

then arguments of the proof of proposition 3.15 still work but instead of individual sub- 
spaces H a = {v{x")T a (x')} for a = (a'; a") = (ai, . . . , a v \ a p+ i, . . . , a r ) one should con- 
sider H a ^ n = {J2\a"\=n V (<x'i<x")( X ") T( Q / ;a //) (x') } . 

Therefore we have 

Proposition 3.17. Proposition 3.15 holds even if /j are not necessary disjoint. 
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3.3.5 Ultrastrong magnetic field case. End 

First of all standard elliptic arguments applied to equation (3.36) yield that under condition 
(0.41) F t ^ h -i T XT(t)u a/ 3 = as r < e' and therefore 

Proposition 3.18. Let conditions (3.16) and (0.36) — (0.37) be fulfilled (with fixed a). Let 
T = /i m with large enough m. Then under condition (0.41) on suppip 

(3.48) \F t ^ T x T (t)^(x)u\ < Cfi~ s Vr : |r| < e'. 

Proposition 3.19. Let conditions (3.16) and (0.36) — (0.37) be fulfilled (with fixed a) . 

(i) Let (f)i be supported in 5(0,1), ^ = 1 in -8(0,2), x be supported in [—1,1], Qk = 

<Pk,M^ 1 T( x " ~ x " 1 £,x~£,' ) w - Then for T == Cs^l log/i| <T <T[ = CqTq and large enough 
constant M estimate (3.7) holds. 

(ii) Let tpk = 4 > k,M^- 1 r(x — x). Then for Tq = f C(fih\ log/i|)2 < T < T x = e/i T and 
large enough constant M and large enough constant M 

(3.49) \F t ^ h -i rX T(t) (1 - faWMx, y, t)<f> x {y) \ < Cfi~ s Vr < c. 

Proof. For a = (3 = a statement (i) follows immediately from (3.36) by standard propa- 
gation method. Then (3.34) implies (i) for all a, (3. Then the similar statement is true for 
u r . 

It implies (3.49) for T^<T <T{; so 0(/iA)U(t)^i = </>(hD t )if) 2 U{T)if)i as T < t < T{. 
In the general case let n = \t/T\. Then 

(3.50) 0(/iA)U(t)^ = U(t - T)^ 2 U(T)^ = . . . 

= 4>(hD t )^ n U(t - nT)Vn-iU(T)V„- 2 • ■ • ^ 2 U(T)^i 

where ipk+i is supported in 2M j u" 1 T- vicinity of suppipk and equal to 1 in in M/i _1 T- vicinity 
of supp^. Then (1 — ip n+2 )(p(hD t )XJ (t)ip\ = which is a statement (ii). □ 

The same arguments lead to 

Proposition 3.20. Let conditions (3.16) and (0.36) — (0.37) be fulfilled (with fixed a). 
Then 

(i) Let microhyperbolicity condition 0.2 which is equivalent now to 

(3.51) £(W S ) > e - C|W a | 
be fulfilled. Let Q 1 = fa^x" - x", £" - t") w and 

(3.52) Q 2 = w^T' 1 ^', x" - x") + (£', x" - x"x" - x")) + T^e'i^ 
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with our standard function w . Then for T G [To, T{] 

(3.53) \F^ h -i T x T {t) (1 - Q2 X )u af3 (x, y, t)Q\J < Cfi~ s Vr < e'. 

(ii) Further, foripi = — x) and ip2 = w(/iT _1 (£, x — x) + T~ 1 e't) withT £ [Tq,Ti]. 

i7ie following estimate holds 

(3.54) |F t ^-i rXT (i)(l -0 3 (z))u(x,j/,t)<Mj/)| < C>~ s Vr < e'. 
Corollary 3.21. Under conditions (3.16), (0.36) — (0.37) and condition 0.2 

(3.55) |iWirXr(*)ru0i| < C^"' Vr : |r| < e' 
/io/rfs /or T £ [To, T]. In particular, 

(3.56) |^Uh-ir(xT(<) - X2b(*))r«0i| < Cfi~ s Vr : |r| < e' . 

Furthermore, for (I, a) >z (2, 0) one can prove easily that singularities propagate along 
trajectories 

(3.57) ^2- = ^J2f ik d Xk V & , V a = X)/ i (2S i + l) + V) Vr:|r|<e' 

k j 

where (f^ k ) = (i 7 ^) -1 but we do not need it. In particular, V a is an integral. For d = 2 
even condition (I, a) >z (2, 0) is not necessary. 

3.4 Propagation. V. Conclusion 

So, under our assumptions to smoothness and microhyperbolicity and choice of e we con- 
sidered cases 

(a) eh~^\ log /i | — ^ < ji < eh~ 1 \ log/i| _1 ; 

(b) eh" 1 ] \ogh\~ 1 < fi < Ch" 1 and #91 = 1; 

(c) eh- 1 ] \ogh\~ 1 <[i< eh" 1 and #9T > 2; 

(d) eh- 1 <fi< Ch- 1 ; 

(e) /j, > eh" 1 and g^ k , Fjk are constant; 

(f) fi > Ch" 1 and condition (0.36)-(0.37) is fulfilled; 

and proved (3. 55), (3. 56), for T = Ce^h] log/i| and Ti = e/i (unless (Z, cr) -< (2,0), and 
#|9t > 2 which is possible in cases (a)-(c) only, when T\ = e/i /_1 | log/i|°"). In cases (a)-(c) 
we knew from section 1 similar result with Tq instead of T and therefore we needed only 
to consider T £ [T ,Tq] here; in cases (d)-(f) we had no help from section 1. 
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4 Short-range propagation and Estimates 
4.1 Goals 

Still, T = Ce _1 /i| \ogh\ is too large. Really, Tauberian method gives us remainder estimate 

(4.1) - sup \F t ^ h -i T (x T {t)rui/;)\ 
while the principal part is 

(4.2) (F t ^ h -, T {x T (t)Tu^))dT. 

We will reduce (4.2) to more explicit form in the next section (where this section approach 
will be crucial as well) but now let us consider (4.1). To decrease it one should increase 
T but generally supremum also grows proportionally T and there is no improvement; the 
following estimate 

(4.3) |F^ h -x T (x T (f)ri^)| < CT(h- 2r + ^h- r ) Vr < c 
rather easily follows from 

(4.4) |0(/iA)(xT(*)r^)| <c(h- 2r + » r h- r y 

However results of the previous section shows that while taking T — T\ one can take xt (t) 
and then (4.1) becomes CT Tf 1 (h~ 2r + fi r h~ r ) . It equals to our dream remainder estimate 

as T = T * = f Ch. Note that even in the smooth case T is larger 21 ). We cannot just take 
T = T * because F t ^ h -i T (xT(t)Tmp) is not negligible for T G [T *,T ]; instead our goal is 
to prove under microhyperbolicity condition that 

(4.5) \F t ^ h -, T (xT{t)Tmf)\ < CT (h~ 2r + ^ h~ r ) x (^)'| log(^)r 

Vr : \t\ < e. 

As (Z, cr) y (1, 1) it would imply that left-hand expression of (4.3) with T = Tq does not 
exceeds its right-hand expression with T = T *. 

21 ) Since "no mollifications needed" has its own pitfalls, we would need actually request Tq = h 1 ^ 5 with 
arbitrarily small exponent 5 > rather than Tq = Ch\ \ogh\. 
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To achieve this goal we apply ^-mollification to At] this will lead to approximation 
error in F t ^f l -i T (j(T(t)Tui/A which we are going to estimate and to (kind of) negligibility 
estimate of it after mollification. Minimizing the sum by rj (which will depend on T) we 
will arrive to (4.5) for non-mollified operator (surely, the original mollification will be still 
here) . 

4.2 Pilot-Model: some classes of pdos 

We start from related simple results (we will need them anyway): 

Proposition 4.1. Let A = A(x,hD) be a (matrix) self-adjoint h-pdo in W with symbol 
a G T l,u and satisfying microhyperbolicity condition on level r = 0: 

(4.6) ((£a)v,v) > e\\v\\ 2 - C\\ (a- t)v\\ 2 \/v V(x,£)- 

Let u be a Schwartz kernel of e~ lh lA . Let ip be supported in the small vicinity ofO. 

(i) Let I > 1. Then for T G [Ch 1 " 5 , T ] with arbitrarily small exponent S > and small 
enough constant To 

(4.7) |^A)xT(0(rw4)l < Chr(±) l -\i + log^r, 

(4.8) |i^_ lrXr (t)(r«^J)| <^ 1 ^(^) i " 1 (l + A)- s |i og ^|- Vr: |r| < e 

with arbitrarily large exponent s. 

(ii) Let I = l,cr > 2. Then for T G [Ch l ~ 5 , To] with arbitrarily small exponent 5 > 
and small enough constant T 

(4.9) \<j>{hD t )xT{t){Tu^ y )\ < Ch~ r \ log (1 + ^)~ s + Ch~ r \ log^r, 

(4.10) \F t ^ h - wX T(t)(Tu4)\ < Ch^ r \ log^| — (1 + + Ch>-r\ bg£|- 

as t : \t\ < e with arbitrarily large exponent s. 

Pilot-model. Before actual proof, let us start with a simple example (actually arising 

asd = 2 22 ): 

22 )ln this case one can apply successive approximations on this stage and to arrive to oscillatory integral 
expression 
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In the case of scalar operator Tr(uip) is something like h d I\ with 
h — J e lh x ^ z ' t q{z) dz and the second expression is h 1 ~ d l2 with 
J 2 = h" 1 JJ x(t/T)e ih ~ lx ( z)t q(z) dzdt. I consider case T < e^h, I > 1. 

Replacing A by its ^-mollification we make an error O in A which leads to error 

0[Th~ li &{rj)) in I 1 . In the same time (multiple) integration by parts with respect to t 
shows that the error in J 2 is does not exceed 

(4.11) Ch~ 2 T 2 J (Th'^X -t\ + l)" 2 dz x #(77) < Ch^T^rj). 

Consider now I\ and I2 with mollified A assuming that \d Zl X\ > e. Multiple integration 
by part with respect to Z\ will transform modified I\ into (If) J e th A ^*<7fc (2) dz where 
qk — 0{rj~ k ~ 1 ${rf)) as k is large enough and therefore modified I± is 

o((J)V*-^(»7))- 

We can treat I2 in the similar way, but we add integration by parts with respect to t as 
well and use of nondegeneracy condition as above. Then modified I2 has the same upper 
estimate. 

Therefore both original I\ and I2 do not not exceed CTh~ 1 ${r)) + (J£\ r]~ k ~ lr d(r]). Picking 

(near) optimal 77 = h/T we get estimates < C(fO ^(^)- ft proves (4.7), (4.8) in this 
special case. 

Proof. 23 ) (I) Note first that due to standard propagation results for T > Che~~ 1 \logh\ 
left-hand expressions in (4.7), (4.8) are negligible. Thus we need to consider only T < 

Ti = f Che~ x \ log/i|. We can consider \ supported in [|, 1] and until the end of the proof 
t, Tfc are positive. 

Consider //-mollification A v of operator A with respect to 24 **: 

(4.12) 1>?7>£. 

. Then the error in the operator (the difference between A v and A) will not exceed Cd^) 
in the sense that 

(4.13) \\(A V - A)\\ < m(r,) 
and identity 

(4.14) e ih ~ HA = e ih ~ HA *> + ik" 1 [ e ift ~ 1 (*-*'> j4 (4 - A v )e ih ~ H ' A ^ dt' 
Jo 

23 ) Of proposition 4.1. 

24 ) lt means that we mollify its Weyl symbol; I remind that original symbol is regular in e-scale. 
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implies that 

(4.15) || (e ih ~ HA - e ih ~ HA ^j \\ < CTh" 1 ^) as \t\ < T, 

(4.16) \XT(t)T((u v - u)i>l) | < Ch'^Td^) 

for standard cut-off ip in £?((x,£), l) C this is exactly what we got in pilot-model. 

On the other hand, microhyperbolicity condition (4.6) and theorem 3.1 of [Brlvr] imply 
that ^(hD^^XT^Tu^y) is negligible provided 

(4.17) h s > rj > ChT~ l \\ogh\ 

(and supp0 C [— e, e]) and picking minimal i] = ChT~ 1 \ \ogh\ we arrive to inequality 

(4.18) \ct>{hD t )T X T{t)u^\ < Ch-*^) 1 " 1 - \logh\ l -° 

due to T > h 1 " 6 . This inequality is almost as good as (4.7): we have an extra factor | \ogh\ l 
in the right-hand expression. 

(II) To get gid of this factor we need more delicate analysis. First, let us iterate (4.14): 

(4.19) e^tA = e ih-HA v y- LiPh'P f e ih-Ht-ti~- -tp)A,r A _ A \ x 

P- J A 

0<p<m-l ^ J 

e ih-HiA v ( A _ A j...( A _ A ^ e ih-H P A V dh . . . dtp + Rm 

with the remainder 

(4.20) R m = —i m h- m [ e ih-Ht-ti—~tm)A v , A _ Aj x 

jh-H^ ^ A _ A ^ ...( A _ ^yh- 1 ^ dti _. dtm _ 

Then the remainder trace norm does not exceed C (Th~ 1 'd(rf)} m h~ d and it is less than the 
right-hand expression of (4.7) for large enoughm (because we gain factor (h/T) l ~ 1 \ \og(h/T)\ 
on each step). 

Now as we take p = 1 we need to consider 

(4.21) e ith ~ lA ^(ih-H{A- A v ))ip. 

Without any loss of the generality one can assume that direction of microhyperbolicity is 
£ = dfa. Let us apply eT-admissible partition with respect to x±: ip' — ^'Ei/^^i)- Here 
and below <p (with different indices) are eT-admissible functions. 
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Due to (4.17) logarithmic uncertainty principle T x rj > C'h\ log/i| holds and therefore 
due to microhyperbolicity (p u e lth A = ip u e ith A (p u (modulo operator, negligible after cut- 
off 4>{hD t )) as |T < t < T where distance between suppy5„ and supp<p„ is exactly of 
magnitude fi~ 1 T. Then to prove (4.7) it is sufficient to estimate properly a trace norm of 
^ v (p u Bil)Lp' l/ x To do this it would be enough to prove that 

(4.22) The operator norm of each "sandwich" (p u Bip'Lp u does not exceed 

Really, then the trace norm of each "sandwich" will be CTh~ r x where T is Xi-size 
of supp<£v, After summation with respect to v factor T disappears and we estimate trace 
norm of (4.21) (after <p(hD t ) cut-off) by /T r M(£) -1 , which implies (4.7), (4.9). 

For p > 2 we replace the integral over A p by the sum of integrals over A p PI {\tj — tj\ < 
e'T,j = 1, . . . ,p} with very small constants e'. We also place partitions elements around 
each copy of the operator exponent, so e lh tjAv is now replaced by tp V2j e lh tjAv (p U2j+1 , 
t = t — ti — ■ ■ ■ — t p . Then either supports of <Pu 2j +i and <^ 2j+2 are disjoint by a distance 
at least ep~ x T for some j and we have a sandwich tPv 2j+1 Bip U2j+2 of the type described in 

(4.22) , or supports of <p Uo and (p U2p+1 are disjoint by a distance at least e/i~ 1 T (and then 
after taking trace we get negligible term) or ip U2j e th tiAn tp v%j+x is negligible operator. 

Since any sandwich has a norm not exceeding Crf \ log r/ \~ a and is accompanied by a 
factor ^, it brings as a result factor (^) 1 | \og{j^)\ l ~ a . 

To estimate the norm of a sandwich of (2.22)-type let us consider the sequence of 
mollifications A p with p = p n = 2 n p , p = max((|),e) and n = l,2,...,n such that 
Pn+i = r\. Then operator B = Y,o<n<n B n witn B n = A Pn - A Pn+1 has a ^(p)^ symbol 

with (p = p n . Then the norm of the sandwich ip u B n (p l/ does not exceed C$(p) x {h/Tp) s with 
arbitrarily large s and summation with respect to n does not exceed the same expression 

h Te - s 

with p = po = h/T; so we get ) (l + — ) which is exactly M s defined in (4.22). 
Thus (4.7), (4.9) are proven. 

(Ill) To prove (4.8) consider first case / > 1. Then plugging (4.19)-(4.20) we see that 
for large enough m we can skip R m . Also we can skip term with p = because it involves 
only ?7-admissible operators. Applying the above arguments we need to estimate 

(4.23) |F^- lrXT (t)(rv/;t)| 

with v p a Schwartz kernel of the term in (4.19) with the same notations as above. Let us 
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consider p — 1 first. One can see easily that 

(4.24) F t ^ h - lT (xAt)e uh - lA ^ = 

(27T)- 1 J XtOK*^ 1(A >- t) dt = 

(27T)- 1 1 XT (t)(-h 2 d 2 t + 7 2 )V^ 1 ^((A, - r) 2 + 7 2 r^t = 

with X(g)(0 = (-<h + l 2 T 2 h~ 2 ) q x(t). Let us pick up 7 = p - 

Let us multiply (4.24) from the left by ip' u which is equal to 1 in eT-vicinity of supp <p„. 
Then modulo negligible we can multiply it from the right by (p v which is equal to in 
eT-vicinity of supp <p„: 

(4.25) F t ^ h -i T (xT(tW„e ith - lA ^ = 

W 1 (fY J mAW* e**- 1 ^ ((A, - r) 2 + 7 2 )"V, 

and therefore 

(4.26) T^- lT ( Xr (t)^e^- 1 ^)Ti^ = 

W'ifY J xwAW^^M^ - ^) 2 + 7 2 )"V, " £^ dt 

sandwich 

where the origin of both factors (p' v is clear: the first one is again due to propagation and 
the second one is just taken of (p v . 

Note that X(q) is a function bounded by C(l + ^) 2q - We already proved that the 
operator norm of "sandwich" does not exceed C^y lp \ log/i|~ pcr (l + s and therefore (4.8) 
for term with p = 1 follows from estimate 

(4-27) Yl W*?M A v ~ r) 2 + 7 2 )"VJ|i < C^V" 2 * 

where ||.||i means trace norm. 

(IV) Let us prove (4.27). Because operator is positive, its left-hand expression does not 
exceed 

(4.28) Tr(((A-r) 2 + 7 2 )"V). 
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Let us consider 7-admissible partition; I remind that 7 > e. Because in .8(21,7) with 
z = (x, £) variation of symbol of A n does not exceed C7, we conclude that (4.29) does not 
exceed 

(4.29) Ch~ r [ Tr'((a(z) -r) 2 + -f 2 )- s dz = C^ r h r [ V (( Wj - r) 2 + 7 2 )" s cfe 

where 0(2;) is a (matrix) symbol of A v and Tr' means a (matrix) trace, uj are eigenvalues 
of a (in order) and integrals are taken over bounded domain inDc M. 27 }, 
Then (4.27) follows from 



(4.30) / {\uj{z) -t\ +~f)- q dz < C7 



l-2g 



as g > I which, in turn, follows from microhyperbolicity. Really, without any loss of the 
generality one can assume that £ = d Zl . Then uniformly with respect to z' = (z 2 , . . . , Zd) 
e < d Zl ujj < C which instantly yields that even if we restrict ourselves by integral over Z\, 
(4.31) still holds. 

So, for term with p = 1 estimate {4.8) is proven. For p > 2 we need to modify our trace 
trick. Consider (p + l)-dimensional integrals representing corresponding terms 

(4.31) (2TiY l h~ p [ XT (t)e ih ~ 1 k t - tl -- tp)A "- tT ) Bx 

J {t>t 1 +-+t p } 

e itih~ l A vB . . . Be it k h-^A V ^ . . . dt 

Replacing 

gife" 1 ((*-*! t p )A n -tr) _ (-tftf* +7 2^^2 + ^yl^hT 1 ((t-ti t p )A^-tr) 

and integrating once by parts with respect to t we gain factor (h/T) 2 ((A? ~~ v) 2 + 7 2 ) 1 
with x, replaced by X(i) — (~d 2 + h~ 2 T 2 j 2 )x, but we also get new terms 



(4.32) ch 1 ^ J xAh + •■• + *„) ((A, - + 7 2 ) 1 e- ih " 1 * r 5e J ' lh " lA "5 • • • x 

Se i*pfc-% dtl . . . dtp 

and may be 

(4.33) ch 2 ~ p J Xt(*i + • • • + - + 7 2 )" 1 e- 4rl(T Be 4tirl ^B ■ 



• x 

„2£ti 1 fa A-T) 
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with some other functions x- These new terms are just lesser-dimensional integrals of the 
same type. As before we frame operator exponents by y?* functions before applying this 
trick. 

(V) Now let / = 1, cr > 1. We need to consider R m and generated by its expression 
which is (m + l)-dimensional integral 



(4.34) (2Try 1 h~ m I XT(t)e lh ~ 1 ( {t - tl -- tm)A ^ tT )Bx 
'{t>ti+-+t m } 

e it lh -^A vB . . . Be it m h-^A dtimmm ^ dt 



We treat it in the same way as before, but in the end we apply estimate 
(4.35) ||(4, + i7rVlli<C7r 



7 1 q as q > 1, 
log 7 1 as q = 1, 



which is proven the same way (but simpler) as (4.27) we see that it does not exceed 
C/i 1 - r (^|log/i|- <T T/i- 1 ) m |log/i| = C7i 1 - r |log/i| m ( 1 - <r > +1 as T > h l ~ S ] for large m it is less 
than the right-hand expression in (4.10). 

Further, one can treat m-dimensional integral (4.33) in the same way. So, estimates 
(4.8), (4.10) are proven completely. □ 

4.3 Case of intermediate magnetic field 

So, assuming that /it* = eh~^ | log h\~^ < /i < [i\ = e/i _1 |log| _1 we consider 

(4.36) T : h 1 - 5 < T < T = C/ih\ \ogh\ 

with arbitrarily small exponent 5 > 0. Note that we need to consider only T > T * = f eo/i" 1 
because for \t\ < T * standard theory takes place. Another restriction from below will 
appear later. 

Proposition 4.2. Let microhyperbolicity condition 0.1 be fulfilled. Let <ft be supported in 
the small vicinity ofO. 

(i) Let I > 1. Then forT G [C/i 1 " 5 ,T ] with arbitrarily small exponent 5 > 

(4.37) MhD t ) X T(t)Fu1>)\ < Ch^^) 1 - 1 ] log^r (1 + A)-, 

(4.38) \F t ^ rX T(t) (Turh) \ < Ch~ d (^) \ log + " s 
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with arbitrarily large exponent s. 

(ii) Let I — 1,<7 > 2. Then for T E [Ch 1 " 6 , Tq] with arbitrarily small exponent 5 > 
estimates 

(4.39) \<f>(hD t ) XT {t)(rui>)\ < 

(4.40) \F t ^ h -t T x T (t)(rwiP)\ < 

W : |r| < e wift arbitrarily large exponent s. 

Proof. As in subsection 4.2, we can consider \ supported in [|, 1] and in what follows t, T^ 
are positive. 

Consider ^-mollification A v of operator At with respect to x", £" with r/ satisfying (4.12). 
Then the error in operator (difference between A v and A) will not exceed Cr] l \ logn\~ s in 
the sense that 

(4.41) \\(Ar,-A)v\\ < CV|logr/|- s (||A;|| + 

and then we just follow the proof of proposition 4.1 with the following modifications: 

(i) we consider /i _1 /i-pdos with respect to x" with "matrix" symbols a; furthermore, 
factor and (fih)~ r comes from the "matrix" trace. 

(ii) Between different copies of B we place not only <p(\) U j but also Q'j where Q'j are 
cut-offs with respect to (x',£') (keeping all in zone {\x'\ + |£'| < C/i -1 } (and also a„ in 
vicinity of r n for all n G 91). 

(iii) We remember that ipr = ip{x' , fi~ 1 D'; x" , fi~ 1 hD") with T l,a symbol and we need 
to mollify ip as well. Then there will be either ipr-q or (ipTrj ~ ^r) factor; in the latter case 
it is possible that it is the only factor containing difference. In the former case we follow 
subsection 4.1, in the latter case we just note that the norm of {ijjT-q — ipr) is 0{rf\ \ogr]\^ a ) 
exactly as in (At v — At) but this time it is not accompanied by rather large factor T/h 
and this makes things much more comfortable. □ 

Then summing over partition of [— T ,T ] \ [T *,T *] we get immediately inequality 

(4.42) \F t ^ h -i T {xT (t) - XT (t))r(m/0| < Ch^ d . 

Really, the left-hand expression here does not exceed the right-hand expression of (4.37) 
integrated over ^ from T * to T which, in turn, does not exceed Ch}~ d \ contribution of 
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extra term (4.39) as I = 1 does not exceed Ch 1 d as well. On the other hand, we know 
from rescaling of the standard results that (for fi < h" 1 ) 

(4.43) \F^ h -i T XT{t)T{ui[>)\ < Ch^ d 

as T = T * and therefore it holds for T = T and it follows from section 1, that (4.43) holds 
for T = T\. I remind that T\ = efi as either #9t = 1 or (I, a) >z (2, 0) and 7\ = efi l ~ l | log fi\ a 
otherwise. 

Then standard Tauberian arguments imply immediately 

Corollary 4.3. Under microhyperbolicity condition 0.1 and fi < h 5 ~ l estimate (1.62) holds. 

In the general case we will need analysis of section 5 to prove this statement. Anyway, 
we need it to get more explicit formula than (1.62) provides. 

4.4 Case of strong magnetic field 

I remind that then \i\ = eh~ l \ log h\~ l < fi < \x\ = Coh~ l and T = Ce~ l h\ \ogh\ where 
either #9T = 1 and 

(4.44) e > C^-^llog/il)*, 
or #91 > 2 and 

(4.45) e > CQT 1 ^ log/i|)5e c ^ |log/l1 

and in the latter case /i < fi\ = €\h~ x . 

In this strong magnetic field case we can apply the same arguments combined with 
propositions 3.7-3.9 and extend proposition 4.2: 

Proposition 4.4. Let eh" 1 ] \ogh\~ 1 < fi < 6ih~ l , T = Ce~ l h\ \ogh\ and either #9T = 1 
and (4-44) hold or #9T > 2 and (4-45) hold. 

Let microhyperbolicity condition 0.1 be fulfilled and let T E [/i 1-5 ,Tb]. Then both state- 
ments (i),(ii) of proposition 4-2 remain true. 

Proof. Just repeating proof of proposition 4.2 without significant difference. We use that 
in operator sense \x'\ + < Cfi" 1 still even if in microlocal sense one should put 
C(fi~ x h\ log h\) 2 instead. □ 

4.5 Case of superstrong magnetic field 

This case is split again into two subcases: fi^ = e\h~ x < fi < fi% = C^h~ l and fi > fi\ = 

c h-\ 
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4.5.1 Case of not-ultrastrong magnetic field 

As before case e\h~ x < \i < ch~ l is just a variation of strong magnetic field case but with 
mandatory #91 = 1 and microhyperbolicity condition 0.2 instead of 0.1. Using the same 
arguments as in the proof of proposition 4.2 we arrive to 

Proposition 4.5. Let eh^ 1 < /i < Ch^ 1 , e = (/i" 1 ^] log/i|)^ and T = Ce~ l h\ log h\ and 
#91=1. 

Lei microhyperbolicity condition 0.2 be fulfilled and let T E [/i 1_<5 ,T ]. T/jen 6ot/j state- 
ments (i),(H) of proposition 4-2 remain true. 

4.5.2 Case of ultrastrong magnetic field 

We need to assume now that either g^ k , Fjk are constant or conditions (0.36)-(0.37) are 
fulfilled. Exactly the same method as before leads us to similar results: 

Proposition 4.6. Let \i > CqH' 1 , e = (ii~ l h\ \ogh\)% , T = Ce' 1 ^ \ogh\ and #91 = 1. Let 

either g^, Fjk are constant and microhyperbolicity condition (0.19) be fulfilled or conditions 
(0.36)-(0.37) and microhyperbolicity condition (0.38) be fulfilled. 

(i) Let I > 1. Then for T E [Cfi 5 h,T ] with arbitrarily small exponent S > 

(4.46) \<P(hD t ) XT (t)(Yu^ y )\ ^^--(^'-'(l + ^-llog^r, 

(4-47) (^-^(^(^^^-(^^(l + ^-liog^- 

as \t\ < e with arbitrarily large exponent s. 

(ii) Let I — 1,(7 > 2. Then for T E \Ch l ~ 5 ,Tq] with arbitrarily small exponent 5 > 
and small enough constant Tq 

(4.48) \<p(hD t ) X T(t)(TmPl)\ < Cf, r h~ r \ log ^|~ CT (l + ^-)~ s + C^h~ r \ log ~|-, 

(4.49) |F^- lrXT (t)(r^t)| < log^| — (1 + + C/i^ 1 -^ log£|- 
as \t\ < e with arbitrarily large exponent s. 

5 Successful approximation and calculations. 

In this section we will replace rather implicit formula (0.48) by a more explicit one and we 
finish the proof of the crucial inequality 

(5.1) \F t ^ h - lT TxT (t)u^\ < Ch l ' 2r + Cfx r h l ~ r . 
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Our tool will be method of successive approximations on rather short interval [—To, To]. 
5.1 Preliminary remarks 

Now even if we have a remainder estimate proven in some cases (see 4.3) but its principal 
part is given by rather implicit formula (0.48) which is transformed into 

(5.2) h- 1 j° (V t _^- lT XT(t) Tr (e^ lu ^r) ) dr 

with T — T\ which can be replaced by T = To (see below) because their difference is 
negligible and ipr = T*tpT. We cannot take T = T * here. 

Surely, microhyperbolicity condition is fulfilled only at levels close to but we do not 
need it otherwise. Really, let us decompose expression (5.2) into sum of 



(5.3) 
and 



(5.4) Mr) (f^-wxtV) Tr(f h ~ HA ^)\ dr 



i n 



where 0, are admissible functions, is supported in [—1, 1] and equal to 1 in [- 
= 9(—t)(1 — 0(r)), L is a small constant. Now we can use microhyperbolicity in the 
first term and replace T = T\ by any T £ [T , T]. On the other hand, one can rewrite the 
second term as 



(5.5) 



h(hDAfat)Tr(e ih - ltA i> T ))) 



t=o 



and then one can replace here T = T\ by any T £ [T *, Ti] where T * = CL~ x h\ log fi\ To. 
One can see easily that (5.5) is of magnitude h~ d + ji r h~ r and therefore we expect the 
principal part of asymptotics of the same magnitude. We assume that microhyperbolicity 
is in the direction dt . 
Note that 



(5.6) 



/ h- l (F t ^ h -, T x T (t)T{u^)) dr = T- l (F t ^ h -i T XiT{t)T(u^) 

J — oo 



r=0 



with Xi{t) = t l x{t)'-> this formula will play very important role below. 
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Remark 5.1. It follows from propositions 4.2, 4.4-4.6 that in their assumptions expression 

(5.6) does not exceed Ch~ d (j;) l \ log(^)| _<J and thus 

(5.7) f h~ l \ (iW lT (x To (t) - Xr(t))r(«V)) dr| < 

+ /*-/!-) (^)'llog(^) I"'. 

In particular for T = T == ee~ 1 h the right-hand expression does not exceed remainder 
estimates we want to prove (in any of our cases). We can even reduce T if smoothness 
allows. In particular, as T = e/i -1 we get remainder estimate Ch~ d (fih) l \ log/i| _CT which is 
a bit better than in theorem 0.3. 

5.2 Case of intermediate magnetic field 

We apply method of successive approximations to calculate both terms (5.3) and (5.4). 
We consider our operator Aq-(x f ,hD';x",hD x rr) as a matrix operator and then we take 
Ar(x' , hD' ,y" , hD x rr) for unperturbed operator A? (simplifying it later). Then one can 
see easily that each next term of successive approximations gets an extra factor Q ± R or 
Q ± R where and Q ± are forward and backward parametrices of hD t — A? and hD t — Aq- 
respectively and by Duhamel principle their operator norms on interval [— T, T] do not 
exceed Th~ l . 

5.2.1 Mollified asymptotics 

Each next term of successive approximations when plugged there adds a factor ||i?|| to its 
estimate from above. 

Really, it follows from the fact that t~ x Q and t~ x Q are operators with norms h~ x 
and factor t either annuls the corresponding term, or is replaced by its commutator with 
(j>L(hD t ) thus releasing factor h. So, with (5.5) any perturbation with norm 0(h s ) would be 
good enough to have only a bounded number of successive approximations to be considered, 
but (xj — Uj) is even better: 

(5.8) Q ± (x j - Vj ) = ( Xj - - g±[A, x^ 

and similarly for and as explained, these parametrices plugged in (5.5) do not increase 
upper estimate, while the norm of [A, Xj] is 0([i~ l h). 

Further, perturbations of the form (At — A) adds factor jjT l \ \ogh\~ a to the correspond- 
ing term in (5.5); we take A as unperturbed operator in this subsubsection. 
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So any term of successive approximation but the first one, leads to term in (5.5) not ex- 
ceeding + \ogh\~ a ) x hr d which, in turn, does not exceed the remainder estimate. 
Therefore 

(5.9) Modulo remainder estimate one can replace in (5.5) operator At by A leading to 
expression 

(5.10) (2 7 r)-V r / i - r J J Mr)d T (rr>(6(T - A(x'\e))Mx'\n)^ dx" de 

where both A(x",£") and ipr(x",£") are considered as "matrices" i.e. as operators in 
auxiliary space H = L 2 (W), 9(t- A(x", £")) is its spectral projector and Tr' is a "matrix" 
trace. 



5.2.2 Non-mollified asymptotics. I 

Unfortunately things are not that good for term in (5.3) because parametrices are "worth" 
of T/h and therefore in view of (5.8) factor (xj — i/j) is "worth" of T 2 h~ 2 x fi~ l h x 
£ _2 yU _1 /i| \ogh\ 2 x fih\\ogh\ 2 (as T = T ) which is not that small, especially for large 

In 2-dimensional case we could always arrange (after appropriate symplectic map 25 - ) ) to 
have (xj —yj) accomplished by extra factor (e + e l ~ 1 \ loge\~ a ) because we considered scalar 
symbol. Here we have essentially matrix symbol and this trick does not work (however it 
will work for ultra-strong magnetic field). A bit larger e does not help much as fi is close to 
hT 1 . Further, perturbation (Aq- — A) is "worth" of Th~ x x fi~ l \ log/i|~ CT x e~ l fi~ l \ log/i| 1_<T 
(again as T = To). In 2-dimensional case we just removed such term at the very beginning. 

However, using methods of section 4 we will be able to "place" factor (^)'| \ogh\~ a into 
estimates of successive approximation terms which will take care of few copies of factor j-. 

So, our goal now is to evaluate properly 

(5.11) F t ^ h -i T XT(t)T (uif) = F t ^ h -i T xr 
as |t| < e. 

The idea of successive approximation method is the following: rewriting down 

(5.12) (hD t - A r )u% = Tih5(t)S(x - y) 
with m ± = u9(±t) and thus 

(5.13) u£ = TihG ± 5(t)5(x - y) 

25 )\Vhich after quantization gives metaplectic transformation. 
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as 

(5.14) (hD t - A T )u% = Tih5(t)5(x - y) + (A T - A r )w ± 
with u ± = ud(±t) and thus 

(5.15) u% = ^hG ± 6(t)5(x -y) + G ± (A T - A r )n ± 
we can iterate the last equation few times: 

(5.16) u% = ^ih (p ± {A r -A T )^g ± 5{t)5{x-y) 

0<fc<m-l 

Tih^iAr-Ar)) G ± 5{t)5{x - y); 

however our latest technique allow usually to take m = 1. Applying ip? from the right we 
get that 

(5.17) u^ = ^ih Yl {0 ± (A T -A r )) k g ± S(t)IC4x,y) 

0<fe<m-l 

Tih^iAr-Ar)) Q^JC^y) 

where K.^(x, y) is the Schwartz kernel of ipr- From now V means that we set both Schwartz 
kernel arguments equal to y' (the same we used in successive approximations) and then 
integrate. 

Proposition 5.2. In frames of proposition 4-2 for T e [h 1 " 6 ,^] 

(5.18) T-'lF^-^XT^T^iAr-A^Sit)^] < 

C^h^(t) l ~ 2 \\oghr. 

Therefore 

(5.19) | f F t ^- lr (xf (t)-Xr*W)r(^ ± (A r -A r )^(t)/C v ,)rfr| < 

J — oo 

c \^ 1 h 1 ~ d as (l,a)h (2,1), 

\fi 1 - l h 1 - d \\ogh\- a as (/, a) -< (2, 1) 

and thus does not exceed remainder estimate of theorem 0. 7. 
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Proof. Writing 

(5.20) A r - A T = ^(xj - y/jRj, Rj = Rj(x', iT x hD'\ x", y", n~ l hD") 

j 

we see that 

(5.21) (g±(A T - A^Q^K^j = 

where we used standard commutation relation (5.8) and fC^j = (xj — Vj)JC^ is the Schwartz 
kernel of [sCj,^]. Note that symbols of A%-\ Rj belong to jF i_1,fT while symbol of [x?,^] 
belongs to fi~ 1 hj rl ~ 1 ' a . 

Therefore due to Duhamel formula the left-hand expression of (5.18) is an absolute value 
of the sum of 

(5.22) - ifi^hT^F^h-irXTit) Tr £ ^R^A^G^t) = 

j 

- i^h-^F^-^XTit) ^Tr x 

3 

( [ jh-Ht-t 1 -t 2 )A TR .jh-H 1 A TA U) e ih-H 2 A T ^ T dti dt \ 

V{i 1 +t 2 <i} ' 

and 

(5.23) T-'F^-^XTit) Tr ^(^^fo^r]) = 

j 

r^^-^XTWVTrf/ e^-^R^-^ixj^^dtX 

For 7] = ChT~ l \ \ogh\ let us make 77-mollification At v of (symbol) of At and also of 
thus we will get also At v - Note that left-hand expression of (5.18) with At, At replaced 
by At v , At v is negligible again due to proposition 3.2. Really, both operators At v , At v 
are microhyperbolic in the same direction. 

Consider the difference between (5.21) for //-mollified and original operators. 

To do this we apply formula (4.19)-(4.20) assuming first that / > 1 and skipping negli- 
gible remainder. Then the difference in (5.22) will be 
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— i/u, 1 h 1 F t _>h-i T XT{t) Tr applied to the sum of the following terms (with constant 
coefficients) 

(5.24) yT 1 h~ 1 ~ k ~ kl ~ k2 f e ^ _1 (*-*i-*2— ■-t k )A TvCie ih-h 1 ATr l (^. . .c^ih-HkArr,^' 

where operators Ck are in some order: one copy of A r or (Ap — Api), one copy of Rj V or 
(Rj V — Rj) and (At — At,) constitute the rest; also ip' = ipr-q or ip' = (ipT — i-'T,) and among 
all the factors there must be at least one factor which is the difference between mollified 
and non-mollified operators. 

Acting as in the proof of proposition 4.1 we can also gain factor 
h 2q T~ 2q ((Ar v — r) 2 + 7 2 ) 9 replacing x by (— d 2 + ^ 2 T 2 h~ 2 ) q x] 2-dimensional integral can 
also appear instead of 3-dimensional (one extra integral comes from Fourier transform). 

Each factor (Aq — At,) or (Aq — At,) is worth of rf\ \ogh\~ u and is accompanied 
by a factor not exceeding T/i" 1 while each factor (AP — A^p) or (Rj V — Rj) is worth of 
f] l ~ 1 \ \ogh\~ a and (ipq — ipT,) is worth of rf \ \ogh\~ a . 

Then using estimate (4.35) with 7 = h/T we derive easily that the difference does not 
exceed CK 1 ^ 1 "^) I logh\ l ~ a and we arrive to estimate different from the required one 
by an extra factor | log h | 1 . 

There are many methods to get rid of this factor; I apply the cheapest one. First, let we 
apply m iterations in successive approximations and note that each extra iteration brings 
factor {i~ x h x T 2 h~ 2 < \ih < | \ogh\~ l and use the same arguments as above; we conclude 
that instead of ^ ± we need to consider (Q ± (Ar — At)) Q ± ; so we got rid of Q ± but paid 
for this by a bit more complicated expression. This expression includes also [xj — Xk,Rj], 
A^} etc but one can see easily that the corresponding terms would satisfy required 
estimate. 

Also, because I > 1 we see from above mollification arguments that all the terms with 
more than one factor of (Aq — Aq- V )^ (A^p — Alp) and (Rj — Rj V ) satisfy required estimate 
and thus we need to consider only terms with exactly one such factor. Further, terms 
containing factor (ipq — ^tv?) satisfy this estimate for sure. 

Further, if we replace in such factor rj by po = h/T we will get terms which can be 
properly estimated and thus in the terms under consideration we can replace (Aq — Ar v ) 
by (A rpo - A Tri ) = E n ®n, [Af - A^) by (1% - Agj = J2 n B' n and (R 3 - R jrj ) by 

(R jpo - R jv ) = J2 n K where B n = (A Tpn - A Tpn+1 ) etc where again p n = 2 n p . 

Now we can just repeat arguments of the proof of proposition 4.1, related to pT x T- 
partition with respect to X\ assuming that microhyperbolicity direction is d^. 

On the other hand, if I — 1, using (4.34) and related integration by part trick (but 
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without /i -1 T-partition), one can see easily that far term in successive approximation satisfy 
estimate in question and thus to get rid of Q ± leaving only Q . 

We can apply mollification and (4.19)-(4.20) skipping terms generated by (4.20); then 
we arrive to estimate which contains in comparison with desired one extra factor |log/i| 
and all the terms save the same as above surely satisfy desired estimate. Now we just apply 
the same /x _1 T-partition with respect to x± etc as above. 

Analysis of (5.23) uses the same technique but is simpler. □ 

5.2.3 Non-mollified asymptotics. II 

Now we are in the shortest-range zone \t\ < T * = h 1 " 6 . 

First of all, for eh 5-1 < ji < eh^\ \ogh\~ l we must estimate F t ^h- 1 rXT^ (uripr) ■ Using 
successive approximation method in its rough form we conclude that modulo term not 
exceeding CTh~ d x fi~ 1 h(T/h) 2 x /i~ 1 /?, 1 ~ d ~ 3<5 we can replace here u by r a\ which is a 
Schwartz kernel of e th t - A . Then calculations show that for T = T * and % — 1 on [— |, |] 

(5.25) F t ^ h - w x T T(u T iJr) = 

(27iy r fi r h 1 - r d T J Tt'(9(t - A(x", D)^r(x", £")) dx"d£" 

modulo negligible term where 6(r — A(x", £")) is the spectral projector of A considered as 
operator in H, ip(x", £") is also operator in M and Tr' is a trace in HI. 

One can see easily that 26 -* for small increment dr eigenvalue luj(x", £") of matrix A(x", £") 
belongs to [r, r + dr] only for (x",£") belonging to the set of measure not exceeding Cdr. 
Since there are x (fih)~ r of such eigenvalues we conclude that expression (5.25) does not 
exceed Ch}~ d . This and the previous arguments imply immediately estimate (4.42) with 
T = T * and thus (4.42) with T = T Q and thus (4.42) with T = 7\. Then standard Tauberian 
arguments imply estimate (1.62). So we have proven 

Proposition 5.3. Under microhyperbolicity condition 0.1 for fi < eh~ 1 \ \ogh\~ 1 estimates 
(442) and (1.62) hold. 

Still we need to get more explicit formula: 

Proposition 5.4. Under microhyperbolicity condition 0.1 for \i < th~ 1 \ log/i| -1 

(5.26) |r(#)(r) - (27r)- r /i r /i- r J Tr'^(r - A T {x", £,"))^ r {x" , £")) dx?d£"\ < 

Cii- l h 1 ~ d [l +e / - 2 |log/i|- <T ). 

26 ^Due to microhyperbolicity condition. 
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Proof. In view of what we proved already we need to prove that 
(5.27) \h~ l f F^,- lr XT(t)r(^ ± (A r - Ar^Sit)^) dr\ 

with T = Tq does not exceed right-hand expression of (5.26). 

Consider expression (5.27) first. Clearly it does not exceed C [i~ 1 h~ d ~ 2 T z and for T = Tq 
it is just C/i _1 /i 1 ~ d ~ 3<5 and therefore the required estimate holds as / < 2; so we need to 
consider case I > 2 only. 

Applying then brute-force successive approximations we conclude that replacing in 

(5.27) Q ± by leads to a much smaller error. Also, note that A? — A? = ~ 

yj)A T (j) + Y,j,k( x o ~ Vj)( x k — yk)R'jk with R'jk e 3 rl ~ 2 ' a . Then using commutator equality 
for all (xj — yj) factors, we conclude that contribution of these terms will be also much 
smaller. In the end we are left with term 

(5.28) -i^J^^h-wXTit) TT^Aru^A^Q^r) 

j 

which is reduced to 

(5.29) i(2n)- r ^-i h -r i*U-i T x T (t) x 

j 

Tr' J (£V, e)A r{j) {x", r)£ v, e)A^\x", v, nv> v, a) 

where we have just matrix- valued symbols a(x",£") = Ar(x",£") and G ± {x",£") and we 
replaced matrix-valued function iJjt(x",£") by scalar-valued ip°(x",£"); this is legitimate 
since ip T (x",£") = if>°(x",£") modulo oj/x" 1 ). 

Using that d^G^x" ,£") = d^g^x", £")A T(j) (x", ^d^G^x", £") and similar identity 
for one can rewrite (5.29) due to trace property and scalarity of ip° as 

(5.30) (27r)- r ^- 1 h- r F t ^ h -t T x T (t)x 

l -{2-nr^- l h- r F t _ h - w x T {t) Tr' J G ± (x", V, O 

with B = — ~ ^r(j) anc ^ ^' = ^(vp- ^ e can rewr ite the first term here as convolution 
with respect to r of Lx(Lt) (with L = Th~ l ) with 

(5.31) (2vr)-> r '- 1 / i - r Tr' d T j 6{t - A r (x", £"))B(x", Vf") dx' 'd? '. 



5 SUCCESSFUL APPROXIMATION AND CALCULATIONS. 



75 



Due to microhyperbolicity and monotonicity of 8 with respect to r this convolution does 
not exceed C\x r ~ l b}~ r x (jih)~ r x (1 + ji 2 ~ l \ log/i| _ °") where the last factor is just upper 
bound of B(x",£") and (fih)~ r comes from trace Tr' and what we got does not exceed the 
right-hand expression of (5.26). 

We can rewrite the second term in (5.31) as convolution with respect to r of Lx(Lt) 
(L = Th~ l ) with 

(5.32) (27r)-V^~ r Tr' J q( t _ A T {x" , dx " d £" 

and this convolution does not exceed the same expression where factor (1 + /i 2- '| log/x| _<T ) 
this time is an upper bound for if}'. □ 



5.2.4 Proof of Theorem 0.5 

Now we can decompose 

(2n)- r fi r h- r j Ti'(e(r - A T (x", ?'))i/>t(x", £")) dx"d^" 

as 

(5.33) (2n)- r f/h- r J Ti'6(t - A (x", £") dx" d^' + E corr 

with 
(5.34) 

E COII d ^ (2n)~ r fi r h' r J Tx'(e(r - A r (x" , £")) - 0{r - A,{x" , £''))) V>V, £") 

(2ny r Sh~ r J Tr>(6(r - A T (x", £")) (M*?, O " W, £"))) dx"di". 

Then making change of variables in integration we can rewrite the first term in (5.33) 
as j £ MW (x,t)i/)(x) dx (using (1.17) and (2.28) for Jacobian). 

Let us represent (5.34) as J £f^ T (x, r)ip(x) dx. While it is clearly the case with the first 
term, in the second term we can rewrite 

(5.35) W^,0-^V.O = 

l<\a\ + \f3\<m 
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with m = \_(l,cr)\; then we can just skip the last term (because the output of it will be 
0(n~ l \ log^h~ d )). Since ^ o ^ = V :i< 7 <| a | +/ 3| Pa^d 1 'ip(x) we can always rewrite the 
second term in the required form as well. What we are lacking is the estimate for these two 
terms. 

Proposition 5.5. Under microhyperbolicity condition 0.1 

(i) E co „ < Cfi- l h- d \ logfi\- a for I < 2; 

(ii) E co „ < Ci2- 2 (i2h) l ' 2 h~ d \ \og(/2h)\- a + C^h 1 -* for2<l< 4; 
(Hi) In particular, E COTT < Cfi~ 1 h 1 ^ d for (I, a) >z (3,0). 

So, / > 2, (I, a) -< (3, 0) is the only case when we cannot skip E COTT . 

Proof. First of all, let us replace A? by A. Due to microhyperbolicity and the fact 
that (Aq — A) is bounded by Cfi~ l \ log/i|~ CT , such operation causes an admissible error 
0(n~ l \ \ogfi\- a h~ d ). 

(i) Replacing in the second term in (5.34) A by Aq we get an error 0(fi~ 2 h~ d ); each 
factor (A — Aq) and ipq — ip° contributes factor /i" 1 into remainder estimate. So we need 
to consider the second term in (5.34) with Aq- replaced by Aq and with (ijjq — replaced 
by ip l } 

(5.36) lT= £ 

\a\ + \j3\=m 

but then after taking trace Tr' we get 0. 

In the first term in (5.34) A = Ao + A\. Consider Aq as an unperturbed operator 
and .4.0 + Ai as a perturbed one. Let us return back to representation (1.62). In view of 
our previous results we should take T = Tq = fih here. Using successive approximation 
procedure and replacing x%(f) by xrif) with T e [h, Tq] we see that m-th term of successive 
approximations gives contribution bounded by Ch~ d (fi~ 1 T/h) m ~ l = C h~ d (T /Tq) 771 " 1 . 

Applying our mollification approach one can prove easily that for T G [h 1 " 5 ,Tq] contri- 
bution of m-th term is bounded by C fi 1 ~ m h~ d (h/T) l+1 ~ m \ \ogh\~ a . Then for m > 1 + 1 the 
sum over t-partition does not exceed Cfi~ l h~ d \ log/z|~ <J . So we can take m = 3 here. On 
interval [— h 1 " 5 , h 1 " 6 ] we just use brute-force estimate of the previous paragraph. 

The first term of successive approximations annihilates with the corresponding expres- 
sion for an unperturbed operator. I remind that Aq = (jj 2 J2j k a jk(j(k + v\ preserves 

spaces H m and Aq = (2 Re/1 2 Y,j,k, P b jk P (j(k( p + J2k h kCk) maps H m into H m+ i © ft m _i 
and therefore after taking trace the second term of successive approximations results in 
for ipr replaced by 
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We need also to consider the similar expressions for ipr replaced by ifj 1 but in this case 
the first term of successive approximations vanishes and the second term has a desired 
estimate. 

(ii) On the other hand, for I > 2, (/, a) -< (3, 0) we need to remember that A = 
Ao + A\ + A2, ipr = ip° + ip 1 + ip 2 and also we need to consider four extra terms (the rest 
definitely is estimated in desired way) which are 27 ) 

(5.37) ih~ l f F t ^ h -i T x T (t)Tv(y2TQ^A 1 g^A 1 g^°)dT, 

(5.38) ih- 1 f F t ^ h - lT x T (t) Tr (V T^A 2 Q^°) dr, 

ih' 1 [ F t ^ h -i T x T (t) Tr (V TGoA^ 1 ) dr, 

J —00 _|_ 

ih- 1 [ F t ^ h - wXT {t) Tr (V TGoi'2) dr 



(5.39) 
(5.40) 



with T = T = fih and Q$ parametrices of hD t — Ao. Here Tr means that we take x = y 
and then integrate; I remind that all operators here depend on y. 

Replacing xrif) by Xrif) and applying the same technique as before we find that such 
modified expressions (5.37)-(5.40) do not exceed C n~ 2 h~ d {h/T) l ~ 2 \ \og(h/T)\~ a . Then af- 
ter summation over t-partition we find that contribution of segment [T, T ] is estimated 
now by C/2~ 2 h^ d (h/f ) l ~ 2 \ \og{h/f ) We take f = e/i" 1 . This is clearly possible as 



5-1 



(5.41) fi < h 

Assuming this, let us consider contribution of [— T, T]. Returning from T reduction to 
initial settings we can rewrite (5.37)-(5.40) (modulo terms, estimated properly) as the 
similar expressions but with Q replaced by parametrices G 1 * 1 of operator A, given by 
(1.65): 

(5.42) A = Y t 9 ik PA + V, 9 ik = 9 ik (v),V = V(y), 



p. = hD . _ V .( y ) _ ^(dkVjMix, - y k ), 



27 ^As we return back to parametrices, cut-off on interval [—To, To] etc. 
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and Aj, i$ are replaced by 



(5.43) 



fi 2 -^b a ( y )P a + ^ H Uy)P a , 




(5.44) 



^2 »- ]a] My)P' 



a:\a\>j 



respectively. 

Note that coefficients of all operators are smooth in e-scale (with e = /i" 1 ) which is 
larger than Ch\ \ogh\ required by the standard theory which is applicable then (after we 
scale x i — ► lix, h \— > Lih, t \—> fit producing that 

- before rescaling contribution of intervals [— e/i _1 , — Ch\ log h\] and [Ch\ log h\, e/i -1 ] are 
negligible and 

- the contribution of interval [— Ch\ log h\,Ch\ log h\] is equal to fi~ 2 X]j>o ^{^Y with 
coefficients Kj which do not depend on fx, h. 

We should not care about terms with j > 1 here because these terms do not exceed 
Cfr l h x ~ d . Therefore, E coxx = K fi~ 2 h~ d + 0(fi~ 2 (fih) l ~ 2 \ log(fih)\~ a + fi~ l ti) (we can alway 
replace it modulo properly estimated). 

Now notice that the "main part" of asymptotics is J £ MW (x)ijj(x) dx which in compar- 
ison with theorem 0.3 for ji x h~? \ log/i| - 5 implies that kq = 0. 

So, under assumption (5.41) statement (ii) of proposition is proven. 

(iii) To get rid of assumption (5.41) let us return to terms (5.37)-(5.40) with x replaced 
by x an d with T > e/i" 1 . 

Each of these expression is equal to fi~ 2 h~ d f(T/h,fih) and therefore the above argu- 
ments imply estimate 



This is really strange inequality because its conditions involve T and h separately. This 
gives us certain flexibility. Namely, let us assume that 



(otherwise everything is fine) and replace h h- > hX, T i— ► TA, \i i— > fiX 1 which affects neither 
[xh nor (h/T). Conditions of (5.45) should be fulfilled now with these modified h,T, li. To 
fulfil the second one we take A = h~ 1 (hT~ 1 )s and this is greater than 1 due to (5.46). 



(5.45) /ih\ log/i| < e, T > h 1 ' 5 , e > C /x 



|/(-,^)| < C(-)^ 2 |log(-)|— . 



(5.46) 



T < h 
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The first of assumptions of (5.45) will survive 28 ) but the third one is now replaced by more 
restrictive e > C fi~ 1 \. To get rid of it let us apply C/i~ 1 A-mollification which leads to the 
error in f(T/h,fih), not exceeding 

cv^a)'- 1 ! logAi-^rc^) 2 + cdi-'xyi log/i-^r^^) 3 . 

This is clearly less than the right-hand expression in (5.45) provided T > h(fih)~ s ' with 
arbitrarily small exponent 5' > (5 depends on it). Now inequality (5.45) is proven under 
this humble assumption. Then one can take T = e/i -1 which clearly satisfies it. 

The rest of arguments of (ii) hold with no alternation. □ 

Remark 5.6. Case I = 2 is missing in the statement of proposition 5.5. 

(i) Using the same arguments as above we find that contribution of zone {T < \t\ < T } 
does not exceed Cfi~ 2 h~ d \ log(yu/i)| 1 ~ cr as a > 1 and Cfi~ 2 h~ d \ log h\ l ~ a as a < 1 and 
C/i- 2 /i- d |log(| \og(fih)\/\ \ogh\) | as a = 1. 

On the other hand, contribution of zone {\t\ < T} will be smaller as a < 1 and 
KQH~ 2 h~ d + 0(/i _2 /i _d | log(/i/i) 1 1 ~ <T ) as cr > 1 and the same arguments as above show that 
k = 0. 

(ii) Further, if either there are no third-order resonances or if g^ k and Fjk are constant, 
then Ai = after proper reduction. I can exploit it properly only for a < 0: one can see 
easily that then E COTI < Cfi~ 2 h~ d \ log/i|~ CT . 

5.3 Case of strong magnetic field 

In the strong magnetic field case one can apply the same arguments in the successive 
approximation method remembering that now e > Cbi^hlloghD* and T < he^ 1 < 
(Hlog/T 1 )^. 

Then propositions 5.2-5.4 and their proofs remain true with this minor modification: 

Proposition 5.7. Under microhyperbolicity condition 0.1 for ji < eh~ l 

(i) Estimates (5.18) and (5.19) hold; 

(ii) Estimates (442), (1.62) and (5.26) hold. 

► Then theorem 0.7 is proven. Also, as = 1 all these results hold for \x < Ch" 1 . 

Furthermore, arguments of proposition 5.5 related to T > T = e/i _1 do not change 
either. However, we need to reconsider contribution of the segment [— T, T] as / > 2. 
Again, we need to consider only the sum of terms (5.37)-(5.40). 

28 )ln fact, it will be fulfilled as long as fih\ \og(h/T)\ < e' . 
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We already know from the proof of proposition 5.5(ii) that under assumptions fih\ log h < 
e, e > Cjj," 1 this sum does not exceed Cfi~ 2 h~ d (fih) l ~ 2 \ \og(fih\)~ a . On the other hand, 
all terms (5.37)-(5.40) with T = T and their sum have the form fi~ 2 h~ d f(fih) and then 
repeating arguments of part (hi) of the proof of proposition 5.5 with A = fi K h K ~ 1 we arrive 
to the same estimate Cfi~ 2 h~ d (fih) l ~ 2 \ \og(jih\)~ a for the sum as fih < e'. Furthermore, 
remark 5.6 remains true as well: 

Proposition 5.8. Under microhyperbolicity condition 0.1 for /i < eh" 1 proposition 5.5 and 
remark 5.6 remain true. 

5.4 Case of superstrong magnetic field 

I remind that in this case we use decomposition (3.17) and proposition 3.15 hold and we 
have equation 



where U = a ,/3ea with #21 < C and furthermore #21 = 1 under condition (0.36)-(0.37) 
and A = (cy4 Q/ g) Q , ]j a g2t having symbol of JF' ,<T . 

I also remind that we replace condition 0.1 by condition 0.2 which is equivalent to 
microhyperbolicity (at level 0) condition to A. Then we can apply successive approximation 
method to this system rather than to the original reduced operator and prove immediately 

Proposition 5.9. Let either eh' 1 < \i < Ch" 1 and = 1 in the microhyperbolicity 
condition 0.2, or (0.35) -(0.19) hold, or (0.36) -(0.38) hold. 
Then 

(i) The following estimates hold: 



(5.47) 



(hD t - A)U = 



U\ t=0 = 5(x' - y') 



(5.48) 




with Tr 







£ 



h and 




(5.50) 




CK-W-0- 2 |lo g/ , 
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forT E [Chn 8 ,Ti], \t\ < e. 
(Hi) Therefore 

(5.51) |r(#)(r) - hT x I h t ^ h -i T ,XT{t)T{i)u)\ dr'\ < C^ x h l ~ r Vr : |r| < e 

J — oo 

for any T E [f ,Tx], |r| < e. 

This immediately imply 
Proposition 5.10. In frames of proposition 5.9 

(5.52) {h- 1 [ (F^ ft _vxr(*)r(^u)) dr'- 

J —oo 

(2n)- r ^ 1 h- r Ti'e(T - A(x",Z")^°(x",Z")dx"dZ"\ < 

Cn r - X h x ~ r + Cfi r h- r e l \ log/i|- CT . 

► Then theorem 0.9 is proven. 
5.5 Condition |V| > eo 

Remark 5.11. (i) Assumption \V\ > eo is automatically fulfilled in the superstrong mag- 
netic field case. In the intermediate and strong magnetic field cases we divided by V in 
microhyperbolicity condition 0.1. However, it is not a problem: we could assume in advance 
that 

(5.53) \V\ + |W| > e 

and then condition 0.1 will be needed in zone {x : |V(x)| > e } only. This approach is not 
applicable to estimates of Eco rr . 

(ii) In the case of a weak magnetic field case we can apply rescaling technique with 
scaling parameter 7 = e\V\. Then h i— > h^~2, ji i— > ^72 as long as 7 > C(/i~ 2 + /i 4 /i 2 ) 
and e = Cfih\\ogh\ 1— > 5 ncw = C^~ l nh\ \ogh^~^ |. Then contribution of B(x,j(x)) to 
semiclassical error (C/x -1 ^ 1 "^) does not exceed C ^L~ l h l ~ d ^ d ~ 4 ^ 2 and multiplying by •y~ d o?7 
and integrating we get C[i~ x h}~ d again. With semiclassical error C fi 1 ~ l h 1 ~ d \ log/i| _fJ we do 
equally good as (Z,cr) >- (2,0). 

For jj, < h~s as \V\ < fi~ 2 we take 7 = Cfi~ 2 and then h 1— > hfi 3 < 1, jui-+ 1 and we are 
in frames of the standard case when assumption \V\ > e was not needed. 
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For < ji < as \V\ < /i 4 /i 2 we take 7 = C/i 4 /i 2 and then h 1— > /i~ 6 /i _1 , t— > // 3 /t, 
and we find ourselves in the intermediate magnetic field case when assumption \V\ > e was 
not needed as well (unless we want to estimate correction terms). 

On the other hand, in the original scale e ncw is replaced by e ncw 7 = e and thus one can 
check easily that approximation error remains C(fih) l \ log h\ l ~ cr h~ d . 
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